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It is well-entrenched folklore that all torsion gravity theories predict observationally negligible 
torsion in the solar system, since torsion (if it exists) couples only to the intrinsic spin of elementary 
particles, not to rotational angular momentum. We argue that this assumption has a logical loophole 
which can and should be tested experimentally, and consider non-standard torsion theories in which 
torsion can be generated by macroscopic rotating objects. In the spirit of action=reaction, if a 
rotating mass like a planet can generate torsion, then a gyroscope would be expected to feel torsion. 
An experiment with a gyroscope (without nuclear spin) such as Gravity Probe B (GPB) can test 
theories where this is the case. 

Using symmetry arguments, we show that to lowest order, any torsion field around a uniformly 
rotating spherical mass is determined by seven dimensionless parameters. These parameters ef- 
fectively generalize the PPN formalism and provide a concrete framework for further testing GR. 
We construct a parametrized Lagrangian that includes both standard torsion-free GR and Hayashi- 
Shirafuji maximal torsion gravity as special cases. We demonstrate that classic solar system tests 
rule out the latter and constrain two observable parameters. We show that Gravity Probe B is 
an ideal experiment for further constraining non-standard torsion theories, and work out the most 
general torsion-induced precession of its gyroscope in terms of our torsion parameters. 

PACS numbers: 04.25.Nx, 04.80.Cc 



I. INTRODUCTION 



(S), defined as follows: 



Einstein's General Theory of Relativity (GR) has 
emerged as the hands down most popular candidate for 
a relativistic theory of gravitation, owing both to its 
elegant structure and to its impressive agreement with 
a host of experimental tests since it was first proposed 
about ninety years ago 1-3]. Yet it remains worthwhile 
to subject GR to further tests whenever possible, since 
these can either build further confidence in the theory or 
uncover new physics. Early efforts in this regard focused 
on weak-field solar system tests, and efforts to test GR 
have since been extended to probe stronger gravitational 
fields involved in binary compact objects, black hole ac- 
cretion and cosmology [3-[3o| . 



A. Generalizing general relativity 

The arguably most beautiful aspect of GR is that it 
geometrizes gravitation, with Minkowski spacetime be- 
ing deformed by the matter (and energy) inside it. As 
illustrated in Figure [T] for the most general manifold 
with a metric g and a connection T, departures from 
Minkowski space are characterized by three geometrical 
entities: non-metricity (Q), curvature (R) and torsion 



"This is the "director's cut" version of the article published in 
Phys. Rev. D November Issue, including extra bonus derivations in 
sections §11, §V and Appendix C. 
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GR is the special case where the non-metricity and tor- 
sion are assumed to vanish identically (Q = S = 0, i.e., 
Riemann spacetime), which determines the connection in 
terms of the metric and leaves the metric as the only dy- 
namical entity. However, as Figure [1] illustrates, this is 
by no means the only possibility, and many alternative 
geometric gravity theories have been discussed in the lit- 
erature [20, 136l - t73l ] corresponding to alternative deform- 
ing geometries where other subsets of (Q, R, S) vanish. 
Embedding GR in a broader parametrized class of the- 
ories allowing non-vanishing torsion and non-metricity, 
and experimentally constraining these parameters would 
provide a natural generalization of the highly success- 
ful parametrized post-Newtonian (PPN) program for GR 
testing, which assumes vanishing torsion [l|-[3| . 

For the purposes of this paper, a particularly inter- 
esting generalization of Riemann spacetime is Riemann- 
Cartan Spacetime (also known as U±), which retains 
Q = but is characterized by non-vanishing torsion. In 
C/4, torsion can be dynamical and consequently play a 
role in gravitation alongside the metric. Note that grav- 
itation theories including torsion retain what are often 
regarded as the most beautiful aspects of General Rela- 
tivity, i.e. general covariance and the idea that "gravity 
is geometry" . Torsion is just as geometrical an entity 
as curvature, and torsion theories can be consistent with 
the Weak Equivalence Principle (WEP). 
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FIG. 1: Classification of spaces (Q,R,S) and the reduction flow. 
Metric-Affine spacetime is a manifold endowed with Lorentzian 
metric and linear affine connection without any restrictions. All 
spaces below it except the Weyl-Cartan space are special cases 
obtained from it by imposing three types of constraints: vanish- 
ing non-metricity tensor Qpivp (Q for short), vanishing Riemann 
curvature tensor Rp V p a (R for short), or vanishing torsion tensor 
Siiv (S for short). A plus sign in a parenthesis indicates a non- 
vanishing quantity from the set (Q, R, S), and a minus sign a van- 
ishing quantity. For example, Riemann spacetime ( 1 ) means 

that Q = S = but R ^ 0. Weyl-Cartan space is a Metric-Affine 
space with vanishing "tracefree nonmetricity" Qp Vp (Q for short), 
defined by Qp, vp = Q^vp — \(trQ)pg vp . The trace of the nonmetric- 
ity is defined by (trQ) M = g vp Q liV p\ thus Q automatically satisfies 
that (trQ)p = (tracefree). Subsets of the classification scheme are 
shown in Fig. 2 of 64], Fig. 1 of [iSfj] and Fig. 5 of Among the 
terms, Einstein- Weyl, Weitzenbock and Minkowski spaces are stan- 
dard, Metric-Affine, Weyl-Cartan, Riemann-Cartan and Riemann 
spaces follow |04l . and we here introduce the terms Weitzenbock- 
Weyl and Minkowski- Weyl space by symmetry. 



B. Why torsion testing is timely 

Experimental searches for torsion have so far been 
rather limited [371 ], in part because most published tor- 
sion theories predict a negligible amount of torsion in the 
solar system. First of all, many torsion Lagrangians im- 
ply that torsion is related to its source via an algebraic 
equation rather than via a differential equation, so that 
(as opposed to curvature), torsion must vanish in vac- 
uum. Second, even within the subset of torsion theories 
where torsion propagates and can exist in vacuum, it is 
usually assumed that it couples onl y to intrinsic s pin, not 
to rotational angular momentum [42l 1 106l Il07| . and is 
therefore negligibly small far from extreme objects such 
as neutron stars. This second assumption also implies 
that even if torsion were present in the solar system, it 
would only affect particles with intrinsic spin (e.g. a gy- 



roscope with net magnetic polarization) jl06l4l08l . Ill3l — 
ll 1 7j j , while having no influenc e on the precession of a 
gyroscope without nuclear spin such as a gyro- 

scope in Gravity Probe B. 

Whether torsion does or does not satisfy these pes- 
simistic assumptions depends on what the Lagrangian is, 
which is of course one of the things that should be tested 
experimentally rather than assumed. Taken at face value, 
the Hayashi-Shirafuji Lagrangian (75j provides an explicit 
counterexample to both assumptions, with even a static 
massive body generating a torsion field — indeed, such a 
strong one that the gravitational forces are due entirely 
to torsion, not to curvature. As another illustrative ex- 
ample, we will develop in Section IIXI a family of tetrad 
theories in Riemann-Cartan space which linearly interpo- 
late between GR and the Hayashi-Shirafuji theory. Al- 
though these particular Lagrangeans come with impor- 
tant caveats to which we return below (see also 126]), 
they show that one cannot dismiss out of hand the possi- 
bility that angular momentum sources non-l ocal torsion 
(see also Table 0. Note that the proof [l06^108t of the 
oft-repeated assertion that a gyroscope without nuclear 
spin cannot feel torsion crucially relies on the assumption 
that orbital angular momentum cannot be the source of 
torsion. This proof is therefore not generally applicable 
in the context of non-standard torsion theories. 

More generally, in the spirit of action=reaction, if a 
(non-rotating or rotating) mass like a planet can generate 
torsion, then a gyroscope without nuclear spin could be 
expected feel torsion, so the question of whether a non- 
standard gravitational Lagrangian causes torsion in the 
solar system is one which can and should be addressed 
experimentally. 

This experimental question is timely because the 
Stanford-led gyroscope satellite experiment, Gravity 
Probe B 1 (GPB), was launched in April 2004 and has 
successfully been taking data. Preliminary GPB results, 
released in April 2007, have confirmed the geodetic pre- 
cession to better than 1%, and the full results, which 
are highly relevant to this paper, are due to be released 
soon. GPB contains a set of four extremely spherical gy- 
roscopes and flies in a circular polar orbit with altitude 
640 kilometers, and we will show that it has the poten- 
tial to severely constrain a broad class of previously al- 
lowed torsion theories. GPB was intended to test the 
GR prediction (76rl8l| that a gyroscope in this orbit pro- 
cesses about 6,614.4 milli-arcseconds per year around its 
orbital angular momentum vector (geodetic precession) 
and about 40.9 milli-arcseconds per year about Earth's 
angular momentum vector (frame-dragging) 2 . Most im- 
pressively, GPB should convincingly observe the frame- 
dragging effect, an arguably still undetected effect of the 



1 http://einstein.stanford.edu/ 

2 These numerical precession rates are taken from the GPB web- 
site. 
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TABLE I: A short list of torsion theories of gravity. The "DOF" in the second column is short for "degrees of freedom". In the 
column Vacuum, "N" refers to non-propagating torsion in the vacuum while "P" means propagating torsion. In the column 
Source, "spin" refers to intrinsic spin while "rotational" means rotational angular momentum. 



off-diagonal metric elements that originate from the ro- 
tation of Earth. Of particular interest to us is that GPB 
can reach a precision of 0.005% for the geodetic preces- 
sion, which as we will see enables precision discrimina- 
tion 3 between GR and a class of torsion theories. 



C. How this paper is organized 

In general, torsion has 24 independent components, 
each being a function of time and position. Fortunately, 
symmetry arguments and a perturbative expansion will 
allow us to greatly simplify the possible form of any tor- 
sion field of Earth, a nearly spherical slowly rotating mas- 
sive object. We will show that the most general possi- 
bility can be elegantly parametrized by merely seven nu- 
merical constants to be constrained experimentally. We 
then derive the effect of torsion on the precession rate of 
a gyroscope in Earth orbit and work out how the anoma- 
lous precession that GPB would register depends on these 
seven parameters. 

The rest of this paper is organized as follows. In Sec- 
tion nn we review the basics of Riemann-Cartan space- 
time. In Section [III1 we derive the results of parametriz- 
ing the torsion field around Earth. In Section I TV! we 
discuss the equation of motion for the precession of a gy- 
roscope and the world- line of its center of mass. We use 
the results to calculate the instantaneous precession rate 
in Section [Vl and then analyze the Fourier moments for 
the particular orbit of GPB in Section fVTl In Section fVIIl 
we show that GPB can constrain two linear combinations 
of the seven torsion parameters, given the constraints on 
the PPN parameters 7 and ai from other solar system 
tests. To make our discussion less abstract, we study 
Hayashi-Shirafuji torsion gravity as an explicit illustra- 



3 GPB also has potential for constraining other GR extensions |82j 
than those we consider in this paper. 



tive example of an alternative gravitational theory that 
can be tested within our framework. In Section [Villi we 
review the basics of Weitzenbock spacetime and Hayashi- 
Shirafuji theory, and then give the torsion-equivalent of 
the linearized Kerr solution. In Section [IXl we generalize 
the Hayashi-Shirafuji theory to a two-parameter family 
of gravity theories, which we will term Einstein-Hayashi- 
Shirafuji (EHS) theories, interpolating between torsion- 
free GR and the Hayashi-Shirafuji maximal torsion the- 
ory. In Section [Xj we apply the precession rate results 
to the EHS theories and discuss the observational con- 
straints that GPB, alongside other solar system tests, will 
be able to place on the parameter space of the family of 
EHS theories. We conclude in Section [XT] Technical de- 
tails of torsion parametrization (i.e. Sect ion [TTT]) are given 
in Appendices [A] & EJ Derivation of solar system tests 
are given in Appendix [Cj We also demonstrate in Ap- 
pendix [D] that current ground-based experimental upper 
bounds on the photon mass do not place more stringent 
constraints on the torsion parameters t\ or t% than GPB 
will. 



After the first version of this paper was submit- 
ted, Flanagan and Rosenthal showed that the Einstein- 
Hayashi-Shirafuji Lagrangian has serious defects [l26j |. 
while leaving open the possibility that there may be 
other viable Lagrangians in the same class (where spin- 
ning objects generate and feel propagating torsion). The 
EHS Lagrangian should therefore not be viewed as a vi- 
able physical model, but as a pedagogical toy model giv- 
ing concrete illustrations of the various effects and con- 
straints that we discuss. 



Throughout this paper, we use natural gravitational 
units where c = G = 1 . Unless we explicitly state other- 
wise, a Greek letter denotes an index running from to 
3 and a Latin letter an index from 1 to 3. We use the 
metric signature convention ( — h H — h). 
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II. RIEMANN-CARTAN SPACETIME 

We review the basics of Riemann-Cartan spacetimc 
only briefly here, and refer the interested reader to Hehl 
et al. 42] for a more comprehensive discussion of space- 
time with torsion. Riemann-Cartan spacetime is a con- 
nected C°° four-dimensional manifold endowed with met- 
ric 57^ of Lorentzian signature and an affine connection 
r M such that the non-metricity defined by Eq. (JTJ) with 
respect to the full connection identically vanishes. In 
other words, the connection in Riemann-Cartan space- 
time may have torsion, but it must still be compatible 
with the metric (g^-x — 0). The covariant derivative of 
a vector is given by 



daV V +Y\ D VP, 



V„V V = d„V v -T p „„V c 



[IV V P 3 



(4) 
(5) 



where the first of the lower indices on r A M(T always corre- 
sponds to the index on V^. 

The full connection has 64 independent components. 
The condition of vanishing non-metricity \7 p g up = 
gives 40 constraints, and the remaining 24 components 
are the degrees of freedom of the torsion tensor. 

In the more familiar case of Riemann spacetime, the 
two conditions = and Q pi/p = imply that the 

connection must be the so-called Levi-Civita connection 
(Christoffel symbol), uniquely determined by the metric 
as 




(6) 



In the more general case when torsion is present, the 
connection must depart from the Levi-Civita connection 
in order to be metric-compatible (^7 p gv P = 0), and this 
departure is (up to a historical minus sign) called the 
contorsion, defined as 



k p = < p \ - r p 



(7) 



Using the fact that the torsion is the part of the connec- 
tion that is antisymmetric in the first two indices (Eq.[3]), 
one readily shows that 

In Riemann-Cartan spacetime, the metric is used to raise 
or lower the indices as usual. 

The curvature tensor is defined as usual, in terms of the 
full connection rather than the Levi-Civita connection: 

= W,* ~ d,T" vX + T" ua T^ x rV^A ■ (9) 

As in Riemann spacetime, one can prove that R p Xvp is a 
tensor by showing that for any vector V p , 

^a v " = \r p x v ,v x - S v fV a V» . (10) 

The Ricci tensor and Ricci scalar are defined by contrac- 
tion the Riemann tensor just as in Riemann spacetime. 



III. PARAMETRIZATION OF THE TORSION 
AND CONNECTION 

The torsion tensor has twenty-four independent com- 
ponents since it is antisymmetric in its first two indices. 
However, its form can be greatly simplified by the fact 
that Earth is well approximated as a uniformly rotating 
spherical object. Throughout this paper, we will there- 
fore Taylor expand all quantities with respect to the di- 
mensionless mass parameter 



m 

3 

r 



(11) 



and the dimensionless angular momentum parameter 



a 

3 

r 



(12) 



where a = J/m is the specific angular momentum , which 
has units of length, and r is the distance of the field 
point from the central gravitating body. Here m and 
J are Earth's mass and rotational angular momentum, 
respectively. Since Earth is slowly rotating (e a 1), 
we will only need to keep track of zeroth and first order 
terms in s a . We will also Taylor expand with respect 
to e m to first order, since we are interested in objects 
with orbital radii vastly exceeding Earth's Schwarzschild 
radius (e m <C l). 4 All calculations will be to first order 
in e m , because to zeroth order in e m , i.e. in Minkowski 
spacetime, there is no torsion. Consequently, we use the 
terms "zeroth order" and "first order" below with respect 
to the expansion in e a . 

We start by studying in section IIII Al the zeroth or- 
der part: the static, spherically and parity symmetric 
case where Earth's rotation is ignored. The first correc- 
tion will be treated in section UlI Bl the stationary and 
spherically axisymmetric contribution caused by Earth's 
rotation. For each case, we start by giving the symmetry 
constraints that apply for any quantity. We then give the 
most general parametrization of torsion and connection 
that is consistent with these symmetries, as derived in the 
appendices. The Kerr-like torsion solution of Hayashi- 
Shirafuji Lagrangian given in Section [Villi is an explicit 
example within this parametrized class. In SectionfVl we 
will apply these results to the precession of a gyroscope 
around Earth. 



A. Zeroth order: the static, spherically and parity 
symmetric case 

This is the order at which Earth's slow rotation is ne- 
glected (e a = 0). For this, three convenient coordinate 



4 These two approximations e m <C 1 and £ a < 1 are highly accu- 
rate for the GPB satellite in an Earth orbit with altitude about 
640 kilometers: e m ~ 6.3 X 10~ 10 and e a =± 5.6 X 10 -7 . 
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systems are often employed - isotropic rectangular co- 
ordinates, isotropic spherical coordinates, and standard 
spherical coordinates. In the following, we will find it 
most convenient to work in isotropic rectangular coordi- 
nates to set up and solve the problem, and then transform 
the result to standard spherical coordinates. 

1. Symmetry Principles 

Tetrad spaces wi th sp herical symmetry have been stud- 
ied by Robertson jll8l | and Hayashi and Shirafuji [75[. 
Our approach in this section essentially follows their 
work. 

Given spherical symmetry, one can naturally find a 
class of isotropic rectangular coordinates (t, x, y, z). Con- 
sider a general quantity 0(x) that may bear upper and 
lower indices. It may or may not be a tensor. In ci- 
ther case, its transformation law 0(x) — ► O'(x') under 
the general coordinate transformation x — * x' should be 
given. By definition, a quantity O is static, spherically 
and parity symmetric if it has the formal functional in- 
variance 

0'(x') = 0(x') 

under the following coordinate transformations (note 
that 0(x') denotes the original function 0(x) evaluated 
at the coordinates x'): 

1. Time translation: t — ► t' = t + to where to is an 
arbitrary constant. 

2. Time reversal: t — > t' = —t. 

3. Continuous rotation and space inversion: 

x^x'eRx, (13) 

where R is any 3x3 constant orthogonal (R'R = I) 
matrix. Note that the parity symmetry allows R 
to be an improper rotation. 



of mass, to keep t\ and t% dimensionless. Below we will 
see that in the context of specific torsion Lagrangians, m 
can be naturally identified with the mass of the object 
generating the torsion, up to a numerical factor close to 
unity. 

We call t\ the "anomalous geodetic torsion" and t% the 
"normal geodetic torsion" , because both will contribute 
to the geodetic spin precession of a gyroscope, the former 
"anomalously" and the latter "regularly" , as will become 
clear in Section [V] and IVII 



3. Torsion and connection in standard spherical coordinates 

In spherical coordinates, the torsion tensor has the fol- 
lowing non-vanishing components: 

V(r)=*x;S ' S re e (r) = S r /(r) = h^, (16) 



where t\ and t% are the same torsion constants as defined 
above. 

The above parametrization of torsion was derived in 
isotropic coordinates, but it is also valid in other spherical 
coordinates as far as the linear perturbation around the 
Minkowski spacetime is concerned. The decomposition 
formula (Eq. [7]), derived from W^g^p — 0, enables one to 
calculate the full connection exactly. However, for that 
purpose the coordinates with a metric must be specified. 
In general, a spherically symmetric coordinate system has 
the line element [83[ 

ds 2 = -h{r)dt 2 + /(r)dr 2 + a(r)r 2 [d(9 2 + sin 2 6>d</> 2 ] . 



2. Parametrization of torsion 

It can be shown (see Appendix A) that, under the 
above conditions, there are only two independent com- 
ponents of the torsion tensor. The non-zero torsion com- 
ponents can be parametrized in isotropic rectangular co- 
ordinates as follows: 

S « = h^x\ (14) 

V = h^&hi-xHji), (15) 

where t\ and ti are dimensionless constants. It is of 
course only the two combinations t\m and t2fn that cor- 
respond to the physical parameters; we have chosen to 
introduce a third redundant quantity m here, with units 



There is freedom to rescale the radius, so-called isotropic 
spherical coordinates corresponding to the choice a(r) = 
f(r). Throughout this paper, we make the common 
choice a(r) — 1, where r can be interpreted as (2tt)~ 1 
times the circumference of a circle. To linear order, 



h{r) = 1 + H-, 



where Ti. and T are dimensionless constants. 

It is straightforward to show that, in the linear regime, 
the most general connection that is static, spherically 
and parity symmetric in Riemann-Cartan spacetime with 
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standard spherical coordinates is as follows: 

H\ m 

ti — 



Y r — 
1 ++ — 



r r = -- — 

TV > 

1 ee 
v r 



Y<t> 



2 y r 2 ' 

7i m 

2^' 

1 2 J r 2 ' 

J 7 m 

2?' 
— r + (J 7 + i 2 )m , 

-r sin 2 8 + (T + t 2 )m sin 2 8 , 
r = 1 

r<j> r ) 

j. 1 m 

4- 

1 0r r t2 r 2 ' 

— sin 6* cos , 

fl = cot 6> . 



(17) 



By "the most general" we mean that any other connec- 
tions are related to the one in Eq. (|17| by the nonlinear 
coordinate transformation law 



dx dx dxi 



T '^ X,) ~ dx« dx'- dx^ Ta ^ ix) + 



dx d 2 x° 



dx a dx' u dx' x ' 
(18) 

Note that the terms independent of metric and torsion 
merely reflect the spherical coordinate system and do not 
represent a deformation of spacetime — in other words, 
the special case t\ = t 2 = H = —T = corresponds 
to the connection for Minkowski spacetime. The case 
ti = t% = and H = —T = —2 corresponds to the 
standard connection for Schwarzschild spacetime in the 
linear regime (r 3> m). 



that the underlying laws of physics are symmetric un- 
der rotations, parity, time translation, and time reversal, 
so that the field configurations for various values of J 
can be related to each other. Specifically, we assume 
that J rotates as a vector, reverses under time-reversal, 
and is invariant under time translation and parity. It 
is then possible to define transformations for the field 
configurations, 0(x\3) — > 0'(x'\J), for these same sym- 
metry operations. Here 0'(x'\3) denotes the transform 
of the field configuration that was specified by J before 
the transformation; O may or may not be a tensor, but 
its transformation properties are assumed to be speci- 
fied. The symmetries of the underlying laws of physics 
then imply that the configurations 0(x\J) are stationary 
and spherically axisymmetric in the sense that the trans- 
formed configuration is identical to the configuration that 
one would compute by transforming J — > J'. That is, 

0'(x'\3)=0(x'\3') 

under the following coordinate transformations: 

1. time translation: t — > t' = t + to where to is an 
arbitrary constant. 

2. Time reversal: —t. 

3. Continuous rotation and space inversion: x — > x' = 
R(x) , i.e. x' is related to x by any proper or im- 
proper rotation. 

Below we will simplify the problem by keeping track 
only of terms linear in J/r 2 = e m e a . 



2. Parametrization of metric 

With these symmetries, it can be shown that the first- 
order contribution to the metric is 



B. First-order: stationary, spherically 
axisymmetric case 

The terms added at this order are due to Earth's ro- 
tation. Roughly speaking, "spherically axisymmetric" 
refers to the property that a system is spherically sym- 
metric except for symmetries broken by an angular mo- 
mentum vector. The rigorous mathematical definition is 
given in Section fill B II Subtleties related to coordinate 
system choices at this order fortunately do not matter in 
the e m -C 1 and e a <C 1 limit that we are interested in. 



1. Symmetry Principles 

Suppose we have a field configuration which depends 
explicitly on the angular momentum J of the central 
spinning body. We can denote the fields generically as 
0{x\3), which is a function of coordinates x and the 
value of the angular momentum vector J. We assume 



gti = gu = -^eijkJ J x ' 



(19) 



in rectangular coordinates x M = (t, x l ), where Q is a con- 
stant, or 

,J 



gu = g<t>t = Q— sin 

r 



(20) 



in spherical coordinates x^ = (i, r, 8, </>) where the polar 
angle 8 is the angle with respect to the rotational angular 
momentum J. The details of the derivation are given in 
Appendix [B] 



3. Parametrization of torsion 

In Appendix [B] we show that, in rectangular coordi- 
nates, the first-order correction to the torsion is 



St 



h 



7 



In spherical coordinates, these first-order torsion terms 
are 



above for the torsion around Earth, we obtain 
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Here /i, • • ■ , /s and twi , . . . , m>b are constants. The latter 
are linear combinations of the former. The details of the 
derivation are given in Appendix [B] We call W\,. . . ,w$ 
the "frame-dragging torsion" , since they will contribute 
the frame-dragging spin precession of a gyroscope as will 
become clear in Section fVl 



o t 



S , 



= tl 



m 

2^2 ' 



q V _ 



O t 

D r<j> 



a <f>_ f m 
°r<l> _r2 2^2' 

ma . 2 
wi— sin 0, 
2r 2 
ma 

W2 -r— sin 9 cos ( 
2r 

ma . 2 
w 3 — sin 9, 
2r z 



(22) 



Sty = w A — 



ma 

u>4 — — sm 9 cos ( 

2r 4 



"'5 



Q 9 

to — 



ma 

2^4 ' 
ma 
W *%5 



cote*. 



All other components vanish. Again, ti, t2, W\,W2, w 3 , 
u>4, w§ are dimensionless constants. 

The calculation of the corresponding connection is 
straightforward by virtue of Eq. ([7]). It is not hard to 
show that, to linear order in a Riemann-Cartan spacetime 
in spherical coordinates, the connection around Earth has 
the following non- vanishing components: 



Around Earth 



We now summarize the results to linear order. We have 
computed the parametrization perturbatively in the di- 
mensionless parameters e m = m/r and e a = a/r. The ze- 
roth order (e a = 0) solution, where Earth's slow rotation 
is ignored, is simply the solution around a static spherical 
body, i.e. the case studied in Section [III Al The first or- 
der correction, due to Earth's rotation, is stationary and 
spherically axisymmetric as derived in Section IIII Bl A 
quantity O to linear order is the sum of these two orders. 
In spherical coordinates, a general line element thus takes 
the form 
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(21) 



where dfl 2 = d9 2 + sin 2 6»d0 2 . Here H, T and Q are 
dimensionless constants. In GR, the Kerr metric [H], Hf| 
at large distance gives the constants Ti = —T = G = —2. 
The result Q = —2 can also be derived more gene rally 
as shown by de Sitter [86| and Lense & Thirring [87f . 
As above, J = ma denotes the magnitude of Earth's 
rotational angular momentum. 

Combining our 0th and 1st order expressions from 
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IV. PRECESSION OF A GYROSCOPE I: 
FUNDAMENTALS 

A. Rotational angular momentum 

There are two ways to covariantly quantify the angu- 
lar momentum of a spinning object, in the literature de- 
noted S M and , respectively. (Despite our overuse of 
the letter 5 1 , they can be distinguished by the number 
of indices.) In the rest frame of the center of mass of a 
gyroscope, the 4- vector 5 M is defined as 



S" = (0,S ) 



(24) 



and the 4-tensor S^ v is defined to be antisymmetric and 
have the components 



S 0i = gu, = 



iO 



where i = x,y,z. So = S x x + S y + S Q z z is the 
rotational angular momentum of a gyroscope observed 
by an observer co-moving with the center of mass of the 
gyroscope. The relation between and S^ v can be 
written in the local (flat) frame as 



5 m = 



'u v S. 



pa > 



(26) 



where i// 1 = dx^/dr is the 4- velocity. 

In curved spacetime, the Levi-Civita symbol is gener- 



alized to e^ pa 



-g where g = detg^. It is 



easy to prove that e^ vpcr is a 4-tensor. Then Eq. 
becomes a covariant relation 



'u v S, 



pa ■ 



(27) 



In addition, the vanishing of temporal components of 
and can be written as covariant conditions as follows: 



0, 
0. 



(28) 
(29) 



In the literature [7y|, Eq. (f2"9"]l is called Pirani's supple- 
mentary condition. 

Note, however, that unlike the flat space case, the 
spatial vectors of 5 M and S pu (denoted by S and 5" 
respectively) do not coincide in the curved spacetime. 
The former is the spatial component of the 4-vector 



5 M , while the latter is historically defined as S 1 = 
e ijk S jk - It follows Eq. (EH) that S and S' differ by 
S = S' [l + 0(m,E/r) + 0{v 2 )\ for a gyroscope moving 
around Earth. 



B. Equation of motion for precession of a 
gyroscope 

To derive the equation of motion for 5 M (or S^ v ) of 
a small extended object that may have either rotational 
angu lar momentum or net spin, Papapetrou's method 
[l20j | should be generalized to Riemann-Cartan space- 
time. This generaliza tion has been studied by Stoeger 
& Yasskin |l06l Il07l | as well as Nomura, Shirafuji & 
Hayashi [l08| . The starting point of this method is the 
Bianchi identity or Noether current in a gravitational the- 
ory whose derivation strongly relies on an assumption of 
what sources torsion. Under the common a ssumption 
that only intrinsic spin sources torsion, both [l06l Il07| 
and [1081 ] drew the conclusion that whereas a particle 
with net intrinsic spin will precess according to the full 
connection, the rotational angular momentum of a gy- 
roscope will not feel the background torsion, i.e. it will 
undergo parallel transport by the Levi-Civita connection 
along the free-falling orbit — the same prediction as in 
GR. 

These results of 106Ml08 | have the simple intuitive in- 
terpretation that if angular momentum is not coupled to 
torsion, then torsion is not coupled to angular momen- 
tum. In other words, for Lagrangians where the angular 
momentum of a rotating object cannot generate a torsion 
field, the torsion field cannot affect the angular momen- 
tum of a rotating object, in the same spirit as Newton's 
dictum "action = reaction" . 

The Hayashi-Shirafuji theory of gravity, which we will 
discuss in detail in Section I Villi raises an objection to 
the common assumption that only intrinsic spin sources 
torsion, in that in this theory even a non-rotating mas- 
sive body can generate torsion in the vacuum nearby [75[ . 
This feature also generically holds for teleparallel theo- 
ries. It has been customary to assume that spinless test 
particles follow metric geodesies (have their momentum 
parallel transported by the Levi-Civita connection), i.e. 
, that spinless particles decouple from the torsion even 
if it is nonzero. For a certain class of Lagrangians, this 
can follow from using the conventio nal variati onal princi- 
ple. However, Kleinert and Pelster [l09l IllOl ] argue that 
the closure failure of parallelograms in the presence of 
torsion adds an additional term to the geodesies which 
causes spinless test particles to follow autoparallel world- 
lines (have their momentum parallel transported by the 
full connection). This scenario thus respects the "action 
= reaction" principle, since a spinless test particle can 
both generate and feel torsion. As a natural extension, 
we explore the possibility that in these theories, a ro- 
tating body also generates torsion through its rotational 
angular momentum, and the torsion in turn affects the 
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motion of spinning objects such as gyroscopes. 

An interesting first-principles derivation of how torsion 
affects a gyroscope in a specific theory might involve gen- 
eralizing the matched asymptotic expansion method of 
|112| |. and match two generalized Kerr-solutions in 
the weak-field limit to obtain the gyroscope equation of 
motion. Since such a calculation would be way beyond 
the scope of the present paper, we will simply limit our 
analysis to exploring some obvious possibilities for laws 
of motion, based on the analogy with spin precession. 

The exact equation of motion for the precession of 
net spin is model dependent, depending on the way the 
matter fields c ouple to the metri c and torsion in the La- 
grangian (see Il2l| ). However, in the 

linear regime that we are interested in here, many of the 
cases reduce to one of the following two equations if there 
is no external non-gravitational force acting on the test 
particle: 



or 



Dr 
Dr 



0, 
0, 



(30) 
(31) 



where D/Dt = (da; M /dr)V M is the covariant differenti- 
ation along the world-line with respect to the full con- 
nection. In other words, the net spin undergoes parallel 
transport by the full connection along its trajectory. 5 

In analog to the precession of spin, we will work out 
the implications of the assumption that the rotational 
angular momentum also precesses by parallel transport 
along the free-fall trajectory using the full connection. 



C. World line of the center of mass 

In GR, test particles move along well-defined trajecto- 
ries - geodesies. In the presence of torsion, things might 
be different. The idea of geodesies originates from two in- 
dependent concepts: autoparallels and extremals 6 . Au- 
toparallels, or affine geodesies, are curves along which the 
velocity vector dx^ /d\ is transported parallel to itself by 
the full connection V . With an affine parameter A, 



the geodesic equation is 



dA 2 



dx^d^ 
dA dA 



(32) 



Extremals, or metric geodesies, are curves of extremal 
spacetime interval with respect to the metric g M „. Since 
ds = [— g llu {x)dx ,Jj dx v ] 1 / 2 does not depend on the full 
connection, the geodesic differential equations derived 



5 If an external non-gravitational force acts on a spinning test par- 
ticle, it will undergo Fermi- Walker transport along its world-line. 
This situation is beyond the interest of a satellite experiment, so 
it will be neglected in the present paper. 

6 This terminology follows Hehl et al. [431 . 



from 5 J ds = state that the 4-vector is parallel trans- 
ported by the Levi-Civita connection. That is, with the 
parameter A properly chosen, 



d 2 xP 
dA 2 " 



fJLlt 



dx p dx u 
dA dA 



= 0. 



(33) 



In Riemann spacetime where torsion identically vanishes, 
Eqs. (|32"f and (|33[) coincide. In a Riemann-Cartan space- 
time, however, these two curves coincide if and only if the 
torsion is totally antisymmetric in all three indices [42| . 
This is because the symmetric part of the full connection 
can be written from Eq. (0) as follows: 



r" 



f 



(liu) 



> V LIU 1 VLLJ 




(34) 



Photons are expected to follow extremal world lines be- 
cause the gauge invariance of the electromagnetic part of 
the Lagrangian, well established by numerous experimen- 
tal upper bounds on the photon mass, prohibits torsion 
from coupling to the electromagnetic field to lowest order 
{42l ]. As a consequence, the classical path of a light ray is 
at least to leading order determined by the metric alone 
as an extremal path, or equivalently as an autoparallel 
curve with respect to the Levi-Civita connection, inde- 
pendent of whether there is torsion. 

On the other hand, the trajectory of a rotating test 
part icle is still an open question in theory. Papapetrou 
[l2dj | claims that, even in GR, a gyroscope will deviate 
from the metric geodesic, albeit slightly. In torsion grav- 
ity theories, the equations of motion for the orbital 4- 
momentum differs more strong ly between different ap- 
proaches [H EhI EqE flllUlT7j , and it is an open ques- 
tion to what extent they are consistent with all classical 
GR tests (deflection of light rays, gravitational rcdshift, 
precession of the perihelion of Mercury, Shapiro time de- 
lay, binary pulsars, etc.). To bracket the uncertainty, we 
will examine the two extreme assumption in turn - that 
world lines are autoparallels and extremals, respectively. 

Only the autoparallel scheme, not the extremal 
scheme, is theoretically consistent, for two reasons. The 
first reason is based on the equivalence of the two ap- 
proaches using the two alternative quantities and 
S'*" to describe the angular momentum. The equiva- 
lence is automatic in GR. In a torsion theory, however, 
Eq. ([50)) and (131[) can be simultaneously valid only if 
the trajectory is autoparallel. This can be seen by tak- 
ing the covariant differentiation of Eq. (T2"T[) . Note that 
De^^/Dr = 0. One finds 



" — >Jp<j — U . 



(35) 



This equation is satisfied if Bu v /Bt = 0, i.e. if the gyro- 
scope world line is autoparallel. If an extremal world line 
is assumed, then one has to make an a priori choice be- 
tween 5 M and S^" , since the precession rates calculated 
using the two quantities will differ. 
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The second reason is that for S^, the condition S^Un 



(Eq. pg])) must be satisfied anywhere along the world 
line. Taking the covariant differentiation for both sides 
of Eq. (|2HI), one finds 



S"Duu/Dt = . 



(36) 



assuming DS^/Dt = 0. Obviously, autoparallels are con- 
sistent with Eq. (|36p . while extremals are not. The same 
argument applies for S^, i.e. taking the covariant differ- 
entiation of both sides of Eq. (|29|) . 

Despite the fact that the extremal scheme is not the- 
oretically consistent in this sense, the inconsistencies are 
numerically small for the linear regime m/r <C 1. They 
are therefore of interest as an approximate phenomeno- 
logical prescription that might at some time in the future 
be incorporated into a consistent theory. We therefore 
include results also for this case below. 



D. Newtonian limit 

In Section Hm we parametrized the metric, torsion and 
connection of Earth, including an arbitrary parameter m 
with units of mass. To give to a physical interpretation, 
the Newtonian limit of a test particle's orbit should be 
evaluated. Obviously, the result depends on whether the 
autoparallel or extremal scheme is assumed. 

In the remainder of this paper, we denote an arbitrary 
parameter with units of mass as Too and the physical 
mass as to. Metric and torsion parameters in accor- 
dance with mo are denoted with a superscript (0), i.e. 

If an autoparallel world line is assumed, using the 
parametrization of equations (|23|) . it can be shown that 
the equation of motion to lowest order becomes 



dv 
df 



,(o) 



H (0) 



(37) 



Therefore Newton's Second Law interprets the mass of 
the central gravitating body to be 



,(0) 



too • (autoparallel scheme) (38) 



However, if — Ti.^ /2 = 0, the autoparallel scheme 
fails totally. 

Similarly, for a theory with extremal world-lines, the 
extremal equation in Newtonian approximation is 



dv 
df 



[-W (0) ] mo 



-e r . 



(39) 



Therefore the physical mass of the body generating the 
gravity field is 



H (0) 



m 



m,Q , (extremal scheme) (40) 



as long as Ti^ ^= 0. For the Schwarzschild metric 
(H (0) = -2), to = to . 

After re-scaling to from too, all metric and tor- 
sion parameters make the inverse re-scaling, e.g. 
t\ = (mo/m) since the combination t\m is the 
physical parameters during parametrization of met- 
ric and torsion. This inverse scaling applies to 
n (o) ^ jr(o) 0(0) _ t (0) ( w (o) w (o) af , welL A natural CQn _ 

sequence of the re-scaling is an identity by definition: 

ti—H/2 = 1, (autoparallel scheme) (41) 
or H = —2 , (extremal scheme) (42) 



V. PRECESSION OF A GYROSCOPE II: 
INSTANTANEOUS RATE 

We now have the tools to calculate the precession of 
a gyroscope. Before proceeding, let us summarize the 
assumptions made so far: 

1. A gyroscope can feel torsion through its rotational 
angular momentum, and the equation of motion is 
either D5 M /Dr = or DS^/Dt = 0. 

2. The world line of a gyroscope is either an autopar- 
allel curve or an extremal curve. 

3. The torsion and connection around Earth are 
parametrized by Eq. and (ggj) . 

With these assumptions, the calculation of the precession 
rate becomes straightforward except for one subtlety de- 
scribed below. 



A. Transformation to the center-of-mass frame 

The precession rate dS/dt derived from a naive ap- 
plication of the equation of motion D5 m /Dt = is the 
rate measured by an observer at rest relative to the cen- 
tral gravitating body. This rate is gauge-dependent and 
unphysical, since it depends on which coordinates the ob- 
server uses; for example, isotropic spherical coordinates 
and standard spherical coordinates yield different pre- 
cession rates. The physical observable is the precession 
rate dSo / dt measured by the observer co- moving with the 
center of mass of the gyroscope, i.e. in the instantaneous 
local inertial frame. 

The methodolog y of transforming S to So was first es- 
tablished by Schiff 76J in which he used the 4-tensor S^ v . 
The basic idea using the 4-vector is as follows. Since 
we are interested in the transformation only to leading 
order in (v/c) 2 and m/r, we are allowed to consider the 
coordinate transformation and the velocity transforma- 
tion separately and add them together in the end. We 
adopt standard spherical coordinates with the line ele- 
ment of Eq. (|21[) . The off-diagonal metric element pro- 
portional to ma/r 2 can be ignored for the purposes of 
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this transformation. Consider a measuring rod in the 
rest frame of the central body. It will be elongated by a 
factor of (l+J 7 m/2r) in the radial direction measured by 
the observer in the center-of-mass frame, but unchanged 
in the tangential direction. The 4-vector transforms 
as dx M ; thus its radial component is enlarged by a factor 
of (1 + Tmjlr) and the tangential components are un- 
changed. This can be compactly written in the following 
form: 



-* — * TTl -t 

So = S + Fy-${S ■ r)r. 



(43) 



Now consider the velocity transformation to the center- 
of-mass frame by boosting the observer along the x-axis, 
say, with velocity v. We have the Lorentz boost from 
gi, = (s°,S x ,Sv,S z ) to Sg = (S °,S x ,So y ,S z ) as 
follows: 



So = l(S° 
S X = l(S x 
S V = S y , 
S Z = S z , 

1 



where 7 = 1/yl — v 2 
S^u^ — gives 



vS x ), 
vS°), 



v 2 /2. 



(44) 
(45) 
(46) 
(47) 

The condition 



S" =v-S = vS x 



which verifies that 5*o = in the center-of-mass frame. 
The spatial components can be written compactly as 



Sq — S — -(S ■ 



v)v . 



(48) 



Combining the coordinate transformation and the veloc- 
ity transformation, we find the following transformation 
from standard spherical coordinates to the center-of-mass 
frame: 



So — S - 



2r 6 



r)r — —{S ■ v)v . 



(49) 



The time derivative of Eq. (|49j) will lead to the expres- 
sion for geodetic precession to leading order , i.e. to or- 
der (m/r)v. To complete the discussion of transforma- 
tions, note that the off-diagonal metric element propor- 
tional to ma/r 2 could add a term of order ma/r 2 to 



Eq. (j4"9"|) , which leads to a precession rate proportional to 
(ma/r 2 )v. Since the leading term of the frame dragging 
effect is of the order ma/r 2 , the leading frame-dragging 
effect is invariant under these transformations, so we are 
allowed to ignore the off-diagonal metric element in the 
transformation. 

The transformation law obtained using the 4-tensor 
guv j g different from using — this is not surprising 
because both descriptions coincide only in the rest frame 
of the gyroscope's center of mass. Schiff [76[ gave the 
transformation law from standard spherical coordinates 
to the center-of-mass frame, using S^: 



So = S' + 



^[S'-(r/r 2 )(r-S')} 



~-[v 2 S'- (v-S')v\. 



(50) 



In taking the time derivative of Eq. (|49|) or (f50|) , one 
encounters terms proportional to dv/dt. Eq. (|37l) or (|39|) 
should be applied, depending on whether autoparallel or 
extremal scheme, respectively, is assumed. 



B. Instantaneous rates 

1. Autoparallel scheme and using 

Now we are now ready to calculate the precession rate. 
In spherical coordinates x M = (t,r,9,<f>), we expand the 
rotational angular momentum vector in an orthonormal 
basis: 

S = S r e r + Sgeg + S^Bcj, . 

In terms of the decomposition coefficients, the 4-vector 
is 

S» = (S° , S 1 , S 2 , S 3 ) = (S°,S r ,S e /r,S 4> /rsm9). 

Applying the equation of motion DS^/Dt = 0, trans- 
forming S to 5*o by Eq. (|49[) and taking the time deriva- 
tive using autoparallels (Eq. [37|) . we obtain the following 
instantaneous gyroscope precession rate: 



d5o 
dt 

where 



= n x sv. 



fl G + Hz 



T U 

V t 2 + 

2 4 



(51) 
(52) 

J.j_(fxtf),(53) 



^i 1 + Mi)(^b ' e r )e r 



-(l + /i 2 )<3j3 



(54) 



Here Iloe = is the angular momentum of Earth, 
where / is Earth's moment of inertia about its poles 
and u>e is its angular velocity. The new effective torsion 
constants are defined so that they represent the torsion- 
induced correction to the GR prediction: 



Ml 

/'2 



(wi 



w 2 - w 3 + 2w 4 + w 5 )/(-3Q) , (55) 
W3+w 6 )/(-g), (56) 



Since t\ — H/2 = 1 in the autoparallel scheme, Eq. ([53 
simplifies to 



n G = + 2t 2 )—^{r* v) 



(57) 



In the literature, the precession due to n G is called 
geodetic precession, and that due to Qp is called frame 
dragging. From Eq. (|53[) . it is seen that geodetic preces- 
sion depends on the mass of Earth and not on whether 
Earth is spinning or not. It is of order mv. The frame- 
dragging effect is a unique effect of Earth's rotation and 
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highlights the importance of the GPB experiment, since 
GPB will be the first to accurately measure the effect of 
the off-diagonal metric element that lacks a counterpart 
in Newtonian gravity. The frame dragging effect is of 
order ma, so it is independent of whether the gyroscope 
is moving or static. In the presence of torsion, we term 
f2<3 the "generalized geodetic precession", and O^ the 
"generalized frame-dragging" . 

2. Extremal scheme and using 

We now repeat the calculation of Section IV B 1} but 
assuming an extremal trajectory (Eq. I39[) when taking 
the time derivative of Eq. (|49|) . obtaining the following 
instantaneous gyroscope precession rate: 

d5Vi ^ 73 tci . 3* _ 
— - = flxSo-h — (S -v f, 58 
at r A 

where O — + Hp . 

n G = ^-W+^^xtO, (59) 

and ftp is the same as in Eq. (15^1) . Since TC = —2 in 
the extremal scheme, Eq. (|59l) is simplified to formally 
coincide with Eq. ([ST]) . 



z(CO E ) 




FIG. 2: A Gravity Probe B gyroscope moves around Earth along 
a circular polar orbit with $o = 7r/2. u)q is its orbital angular 
velocity and u)e is Earth's rotational angular velocity around the 
z-axis. 



Extremal scheme and using 



In spherical coordinates, S^ v satisfies 



S — —Sj,, s 
r v 



23 



1 



~ S T 5 S 



31 



-Sg, (60) 



where S r , S g , are the components of S 
cal coordinates, i.e. S — S e r + S eg + S\i 



in sphcri- 
. We now 

repeat the calculation of Section IV B II assuming an ex- 
tremal trajectory (Eq. l39|) and the S^-based precession 
of Eq. ([50]) when taking the time derivative of Eq. (|49|) . 
obtaining the following instantaneous gyroscope preces- 
sion rate: 

^ = nxSo + h^rx (vxS ), (61) 

where 51 = flc + Of • 

D.G and (If are the same as in equations (f59| and (|54[) . 
respectively. 

In both cases using extremals, the precession rates 
have anomalous terms proportional to ti; see Eq. ([58]) ) 
and We call these terms the "anomalous geode- 

tic precession". These anomalies change the angular 
precession rate of a gyroscope, since their contributions 
to dS'o/dt are not perpendicular to So- This is a phe- 
nomenon that GR does not predict. Meanwhile, £2 con- 
tributes to modify only the magnitude and not the direc- 
tion of We therefore term t\ the anomalous geodetic 



torsion and ti the normal geodetic torsion. The torsion 
functions wi,. . . ,u>5 contribute to the generalized frame- 
dragging effect via the two combinations fix and /12 , and 
we therefore term them "frame-dragging torsions" . 



4- Autoparallel scheme and using S^ v 



Repeating the calculation of Section IV B II using the 
S^-based precession rule of Eq. ([50|) gives the exact same 
instantaneous precession rate as in Section [VB 11 This is 
expected since these two precession rules are equivalent 
in the autoparallel scheme. 



VI. PRECESSION OF A GYROSCOPE III: 
MOMENT ANALYSIS 



GPB measures the rotational angular momentum So 
of gyroscopes and therefore the precession rate dSo/dt 
essentially continuously. This provides a wealth of in- 
formation and deserves careful data analysis. Here we 
develop a simple but sensitive analysis method based on 
Fourier transforms. 
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Fourier transforms 



The angular precession rate is 



The Gravity Probe B satellite has a circular polar orbit 
to good approximation 7 , i.e. the inclination angle of the 
orbital angular velocity ojo with respect to the Earth's 
rotation axis (z-axis) is do — tt/2. Hence the orbital 
plane is perpendicular to the equatorial plane. Let the 
y-axis point along the vector Coo and let the x-axis be 
perpendicular to the y-axis in the equatorial plane so that 
the three axes {x,y,z} form a right-handed coordinate 
basis as illustrated in Figured] A gyroscope at a point P 
is marked by the monotonically increasing angle tp with 
respect to z axis. The polar angle of the point P can be 
regarded as a periodic function of tp: 



3m 



(68) 



where wo = ^>oV is the orbital angular velocity and uje — 
loez is the rotational angular velocity of Earth. Here the 
"biases" relative to the GR prediction are defined by 



W = -(1+^+2*2 + 1^), 
K = ^(l + 3Mi-2/i 2 ), 



(69) 
(70) 



-G) 



[1 + (wi +w 2 -w 3 - + w 5 )/Q] , 



6(tp) 



, < <p < 7T 
if , 7T < if < 2lT 



(62) 



where the constant rj reflects the different assumptions 
that we have explored, and takes the following values: 



So for a particular circular polar orbit, dSa/dt(f, if) can 
be regarded as a periodic function of tp, where ro is 
the fixed radius, allowing us to write dSo/dt(f, v) = 
dS /dt(tp). 

Now define the Fourier moments of the precession rate 

as 



i nds 

a « = ^~ / -jr m cos ntpdtp , 
2tt J dt 

- (tp) sin ntpdtp , 



1 [ 2 * dS 



(64) 
(65) 



2tt J dt 



where n = 1, 2, . . ., so that we can write 



dS 

— — (tp) = 3q + 2 (a n cos ntp + b n sin ntp) . (66) 

n—1 



B. Average precession 

We now write equations ([EE]), (EH), dM]), (|57)l . 
and (|6ip explicitly in terms of tp and perform the Fourier 
transforms. The average precession in the three calcula- 
tion schemes above can be compactly written as follows: 



using autoparallels 
-1 using and extremals (71) 
-1 using and extremals 

From the above formulas, we see that the three schemes 
give identical results when t\ = 0. 

For comparison, GR predicts the average precession 
rate 



where Q c ff 



1. 



-£jr(<P) ) = ^eff X S , 



3m 
2r~t 



(72) 



i.e. , u t = Uft 

It is important to note that torsion contributes to the 
average precession above only via magnitudes of the pre- 
cession rates, leaving the precession axes intact. The 
geodetic torsion parameters t\ and t 2 are degenerate, 
entering only in the linear combination corresponding 
to the bias bt- The frame-dragging torsion parameters 
Wi, . . . , ii>5 are similarly degenerate, entering only in the 
linear combination corresponding to the bias 6^. If for 
technical reasons, the average precession rate is the only 
quantity that GPB can measure, then only these biases 
can be constrained. 



C. Higher moments 



a 



dS 
dt 



(tp) ) = (l cB x So 



(57) Interestingly, all higher Fourier moments vanish except 

for n = 2: 



7 The actual GPB orbit has an orbital eccentricity of 0.0014 and 
an inclination of 90.007° according to the Fact Sheet on the 
GPB website. These deviations from the ideal orbit should cause 
negligible (< 10 -5 ) relative errors in our estimates above. 



-SQIuJe ^ x „ z , m , 

a 2 = — ^-g — (1 + \i\)z x i>o + Vti— wowo z + b x 



4r 



i* —3QIue m n » , j- m lex- a z - 

b 2 = — ~ q — (l + ni)x x bo + rjt-L- — loo(oo x - b z 
8rg 4r 
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Here we use the notation So 1 = Sq ■ i, where i denotes 
the x, y and z axes. 

For comparison, GR predicts the following second mo- 
ments (moments with m = 1 and m > 2 vanish): 



n 2 



3iW „ ^ 



ZIuje 
4rg 



(74) 
(75) 



Technically, it may be difficult to measure these second 
moments because of the extremely small precession rate 
per orbit. However, if they could be measured, they 
could break the degeneracy between t\ and t 2 : could 
be measured through the anomalous n — 2 precession 
moment (the second term in Eq. I|73p). The sign am- 
biguity of t\ is due to the relative sign difference be- 
tween the two schemes using extremals and S^ u versus 
S^ . The degeneracy between wi , . . . , could be allevi- 
ated as well, since the linear combination /Ui (defined in 
Eq. ([55])) could be measured through the correction to 
the normal n = 2 precession moment (the first term in 
Eq. (jT3"|) ). By "anomalous" or "normal", we mean the 
term whose precession axis has not been or already been, 
respectively, predicted by GR. In addition, the anoma- 
lous second-moment terms cannot be expressed as the 
cross product of So and an angular velocity vector. 



VII. CONSTRAINING TORSION 
PARAMETERS WITH GRAVITY PROBE B 

The parametrized Post-Newtonian (PPN) formalism 
has over the past decades demonstrated its success as 
a theoretical framework of testing GR, by embedding 
GR in a broader parametrized class of metric theories 
of gravitation. This idea can be naturally generalized 
by introducing more general departures from GR, e.g. 
torsion. For solar system tests, the seven torsion pa- 
rameters derived in Section IIIII define the torsion exten- 
sion of the PPN parameters, forming a complete set that 
parametrizes all observable signatures of torsion to lowest 
order. 

However, most of existing solar system tests cannot 
constrain the torsion degrees of freedom. Photons are 
usually assumed to decouple from the torsion to pre- 
serve gauge invariance (we return below to the experi- 
mental basis of this), in which case tests using electro- 
magnetic signals (e.g. Shapiro time delay and the de- 
flection of light) can only constrain the metric, i.e. the 
PPN parameter 7, as we explicitly calculate in Appendix 
IC II and Appendix IC 21 Naively, one might expect that 
Mercury's perihelion shift could constrain torsion param- 
eters if Mercury's orbit is an autoparallel curve, but 
calculations in Appendix IC 41 and Appendix IC 51 show 
that to lowest order, the perihelion shift is nonetheless 
only sensitive to the metric. Moreover, PPN calcula- 
tions [|[ show that a complete account of the perihelion 



shift must involve second-order parameters in m/r [e.g. 
the PPN parameter 0) , which are beyond our first-order 
parametrization, as well as the first-order ones. We there- 
fore neglect the constraining power of Mercury's perihe- 
lion shift here. In contrast, the results in Section lYl Bl 
show that Gravity Probe B will be very sensitive to tor- 
sion parameters even if only the average precession rates 
can be measured. 

We may also constrain torsion with experimental upper 
bounds on the photon mass, since the "natural" exten- 
sion of Maxwell Lagrangian (9 M — > V M using the full con- 
nection) breaks gauge invariance and introduces anoma- 
lous electromagnetic forces and a quadratic term in 
that may be identified with the photon mass. In Ap- 
pendix [D] we estimate the constraints on the torsion pa- 
rameters t\ and t 2 from the measured photon mass limits, 
and show that these ground-based experiments can con- 
strain t\ or t 2 only to a level of the order unity, i.e. , not 
enough to be relevant to this paper. 

In Appendix [Cl we confront solar system tests with the 
predictions from GR generalized with our torsion param- 
eters. In general, it is natural to assume that all metric 
parameters take the same form as in PPN formalism 8 , 
i.e. 1 



H 
T 



-2. 



= 2 7: 



.(l+ 7 + -ai). 



(76) 
(77) 

(78) 



Therefore, Shapiro time delay and the deflection of 
light share the same multiplicative bias factor [T — 
H)/4 = (1 + 7 )/2 relative to the GR prediction. The 
analogous bias for gravitational redshift is unity since 
(Az//z/)/(A^/^)( Gfl ) = -H/2 = 1. In contrast, both the 
geodetic precession and the frame-dragging effect have a 
non-trivial multiplicative bias in Eqs. (|6"9")) and (|7D|) : 



bt = 



b u = 



2 7 ) 



\(2t 2 



1 , 1.1. 

-(1 + 7+ ~ai) - - (wi + w 2 - w 3 

~2W4 + w 5 ) . 



(79) 



(80) 



We list the observational constraints that solar system 
tests can place on the PPN and torsion parameters in 
Table |TT] and plot the constraints in the degenerate pa- 
rameter spaces in Figure [3] We see that GPB will ulti- 
mately constrain the linear combination t 2 + -Tft\ (with 77 
depending on the parallel transport scheme) at the 10 -4 
level and the combination ui\ +u>2 — W3 — 2w4 + w§ at the 
1% level. The unpublished preliminary results of GPB 
have confirmed the geodetic precession to less than 1% 



8 This may not be completely true in some particular theories, e.g. 
H. 7^ —2 in Einstein-Hayashi-Shirafuji theories in the autoparallel 
scheme, shown in Table ITVl 
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level. This imposes a constraint on fa + -j- 1±\ < 0.01. 
The combination w± + W2 — W3 — 2w^ + w$ cannot be 
constrained by frame-dragging until GPB will manage to 
improve the accuracy to the target level of less than 1 
milli-arcsecond. 



VIII. LINEARIZED KERR SOLUTION WITH 
TORSION IN WEITZENBOCK SPACETIME 

So far, we have used only symmetry principles to de- 
rive the most general torsion possible around Earth to 
lowest order. We now turn to the separate question of 
whether there is any gravitational Lagrangian that actu- 
ally produces torsion around Earth. We will show that 
the answer is yes by exploring the specific example of the 
Hayashi-Shirafuji Lagrangian [75| in Weitzenbock space- 
time, showing that it populates a certain subset of the 
torsion degrees of freedom that we parametrized above 
and that this torsion mimics the Kerr metric to lowest 
order even though the Riemann curvature of spacetime 
vanishes. We begin with a brief review of Weitzenbock 
spacetime and the Hayashi-Shirafuji Lagrangian, then 
give the linearized solution in terms of the seven param- 
eters t±, £2, wi, . . . , W5 from above. The solution we will 
derive is a particular special case of what the symme- 
try principles allow, and is for the particularly simple 
case where the Riemann curvature vanishes (Weitzenbock 
spacetime). Later in Section ITXl we will give a more gen- 
eral Lagrangian producing both torsion and curvature, 
effectively interpolating between the Weitzenbock case 
below and standard GR. 

We adopt the convention only here in Section lVIIII and 
Section IIXI that Latin letters are indices for the internal 
basis, whereas Greek letters are spacetime indices, both 
running from to 3. 

A. Weitzenbock spacetime 

We give a compact review of Weitzenbock space- 
time and Hayashi-Shirafuji Lagrangian here and in Sec- 
tion IVIII Bl respectively. We refer the interested reader 
to their original papers [74|, [75j for a complete survey of 
these subjects. 

Weitzenbock spacetime is a Riemann- Cart an space- 
time in which the Riemann curvature tensor, defined in 
Eq. vanishes identically: 

#W(r) = o. (si) 

Figure Q] illustrates how Weitzenbock spacetime is related 
to other spacetimes. 

Consider a local coordinate neighborhood of a point 
p in a Weitzenbock manifold with local coordinates x^ . 
Introduce the coordinate basis = {(d/dx 11 ) } and 

the dual basis = {{dx^) p }. A vector V at p can be 

written as V — V^E ll . The manifold is equipped with 



an inner product; the metric is the inner product of the 
coordinate basis vectors, 

g(E,j,,E v ) = g{E v ,E^) = . 

There exists a quadruplet of orthonormal vector fields 
efe(p), where ek(p) = e fe ^(p)£ M , such that 

g(e k ,ei) = g liV e^ L e l v = m , (82) 

where rjki = diag(— 1, 1, 1, 1). There also exists a dual 
quadruplet of orthonormal vector fields e k (p), where 
e k (p) — e k (p)E^, such that 

e k "e\=S^, e k »e\ = 8 k l . (83) 
This implies that 

Vkie k fl e l v = g„ v . (84) 

which is often phrased as the 4x4 matrix e (a.k.a. the 
tetrad or vierbein) being "the square root of the metric" . 

An alternative definition of Weitzenbock spacetime 
that is equivalent to that of Eq. (fcTTj) is the requirement 
that the Riemann-Cartan spacetime admit a quadruplet 
of linearly independent parallel vector fields e fe M , defined 
by 9 

= d ll e k v +T\ x e k x =0. (85) 
Solving this equation, one finds that 

F\u = e k %e k v , (86) 
and that the torsion tensor 

Sj = le k \d,e\-d„e k J. (87) 

This property of allowing globally parallel basis vector 
fields was termed "teleparallelism" by Einstein, since it 
allows unambiguous parallel transport, and formed the 
foundation of the t orsion theory he termed "new general 
relativity" |Mo|. 

A few additional comments are in order: 

1. It is easy to verify that the first definition of 
Weitzenbock spacetime (as curvature-free, i.e. via 
Eq. (|8 1 [) ) follows from the second definition — one 
simply uses the the explicit expression for the con- 
nection (Eq. 156")) . It is also straightforward to verify 
that = using Eq. ([84]) and (|85|) . 

2. Eq. (|86[) is form invariant under general (space- 
time) coordinate transformations due to the non- 
linear transformation law (Eq. (JTSj) ) of the con- 
nection, provided that e k and e k fl transform as 
a contravariant vector and a covariant vector, re- 
spectively. 



9 Note that Hayashi and Shirafuji 75] adopted a convention where 
the order of the lower index placement in the connection is op- 
posite to that in Eq. 1851 . 
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Effects 


Torsion Biases 


Observ. Constraints 


Remarks 


Shapiro time delay 


At/At( GR ) = (l + 7)/2 


7 - 1 = (2.1 ± 2.3) X 1CT 5 


Cassini tracking 122] 


Deflection of light 


5 /g(GR) = (1 + 7 )/2 


7 - 1 = (-1.7 ±4.5) X 10" 4 


VLBI 123] 


Gravitational redshift 


(Ais/iy)/(Au/v) < - GR 1 = 1 


no constraints 




Geodetic Precession 


n G /n<g R) = b t 


(7 - 1) + (t 2 + Jfiti) < 1.1 x 10~ 4 


Gravity Probe B 


Frame-dragging 


u F /af H) =6 M 


|(7 — 1 + j«i) - ^(wi + u>2 - wz - 2««4 + ws)\ < 0.024 


Gravity Probe B 



TABLE II: Constraints of PPN and torsion parameters with solar system tests. The observational constraints on PPN param- 
eters are taken from Table 4 of 01- Unpublished preliminary results of Gravity Probe B have confirmed geodetic precession 
to better than 1%, giving a constraint K7 — 1) + (t2 + ^pti)| i$ 0.01. The full GPB results are yet to be released, so whether 
the frame dragging will agree with the GR prediction is not currently known. The last two rows show the limits that would 
correspond to a GPB result consistent with GR, assuming an angle accuracy of 0.5 milli-arcseconds. 




FIG. 3: constraints on the PPN parameters (7, qi) and torsion parameters (ti, £2, mi • • • ^5) from solar system tests. General 
Relativity corresponds to the black dot (7 — 1 = ai = all torsion parameters = 0). Left panel: the shaded regions in the 
parameter space have already been ruled out by the deflection of light (orange/grey) and Shapiro time delay (yellow/light grey). 
Gyroscope experiments are sensitive to torsion parameters. If the geodetic precession measured by Gravity Probe B is consistent 
with GR, this will rule out everything outside the hatched region, implying that —1.5 x 10 -4 < £2 + < 1.1 x 10 -4 (assuming 
a target angle accuracy of 0.5 milli-arcseconds). The unpublished preliminary results of Gravity Probe B have confirmed the 
geodetic precession to better than 1%, giving a constraint [ia + 4f ii| ^ 0.01. Right panel: the shaded regions in the parameter 
space have already been ruled out by Shapiro time delay combined with lunar laser ranging experiment (yellow/light grey). 
Lunar laser ranging constrains |ai| < 10 -4 ; 3J. If the frame-dragging effect measured by Gravity Probe B is consistent with 
GR, this will rule out everything outside the hatched region, implying that + W2 — 11)3 — 2u>4 + w 5 \ < 4.8 x 10~ 2 . 



3. The Weitzenbock spacetime preserves its geom- 
etry under global proper orthochronous Lorentz 
transformations, i.e. a new equivalent quadru- 
plet of parallel vector fields e' is obtained by a 
global proper orthochronous Lorentz transforma- 
tion, eV =A' fcei ". 



B. Hayashi-Shirafuji Lagrangian 

The Hayashi-Shirafuji Lagrangian |75| is a gravita- 
tional Lagrangian density constructed in the geometry 
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of Weitzenbock spacetime 10 . It is a Poincare gauge the- 
ory in that the parallel vector fields e k (rather than the 
metric or torsion) are the basic entities with respect to 
which the action is varied to obtain the gravitational field 
equations. 

First, note that the torsion tensor in Eq. (|57|) is re- 
ducible under the group of global Lorentz transformation. 
It can be decompos ed in to three irreducible parts under 
this Lorentz group [119I ] 11 . i.e. into parts which do not 
mix under a global Lorentz transformation: 



1 



1 



a A 

D p\ > 



(88) 
(89) 
(90) 



Here e^ pa = y/^e„ vpa and e^ CT = e^/^/y are 4- 
tensors, and the Levi-Civita symbol is normalized such 
that £oi23 = — 1 and e 0123 = +1. The tensor tx^v satisfies 
t\nu = tfiXu, g^txpv — g X ^t\pv = 0, and tx^u + t^x + 
t^xp — 0. Conversely, the torsion can be written in terms 
of its irreducible parts as 



the connection) is not an independent variable as in some 
standard torsion theories [42| . Instead, the torsion is ex- 
clusively determined by the tetrad via Eq. (|87|) . The 
resultant field equation is 

^-G^({ }) + V A F^ A + v x F^ x + - \g pv L 2 = . 
2k 2 

(93) 

Here the first term denotes the Einstein tensor calculated 
using the Levi-Civita connection, but the field equation 
receives important non-Riemannian contributions from 
torsion through the other terms. The other tensors in 
Eq. (|93|) are defined as follows: 



F 



(tf u ' x - t» Xv ) + c 2 (g^v x - g pX v u ) 



H n» = 2S^ p F pa u - S ap "F t " pa , 



(94) 

(95) 
(96) 



Since torsion is the first derivative of the tetrad as per 
Eq. (f87|) . the field equation is a nonlinear second-order 
differential equation of the tetrad. Consequently, the 
tetrad (hence the torsion) can propagate in the vacuum. 



S V pX — 7^{txp,v — t\vti) + -^{g\pVu 



gxuVfj.) + exuvpC 



(91) 

In order that the field equation be a second-order dif- 
ferential equation in (so that torsion can propagate), 
the Lagrangian is required to be quadratic in the tor- 
sion tensor. In addition, the Lagrangian should be in- 
variant under the group of general coordinate transfor- 
mations, under the global proper orthochronous Lorentz 
group, and under parity reversal in the internal basis 
(e — ► eo, e a — > — e a ). Hayashi and Shirafuji suggested 
the gravitational action of the following form 75] : 



Ig = J dW=s[^R({}) + cit 
+c 2 v^v^ + c 3 a M a M ] , 



(92) 



where ci,C2,C3 are three free parameters, R({}) is the 
scalar curvature calculated using the Levi-Civita con- 
nection and k = 87rG/c 4 . The vacuum field equations 
are obtained by varying this action with respect to the 
tetrad e v and then multiplying by rj e^- Note that 
in Hayashi-Shirafuji theory, the torsion (or equivalently, 



The Hayashi-Shirafuji t heory differs from the teleparallel gravity 
theory decribed in |l04t jl05j , which is argued to be fully equiv- 
alent to GR. 

Note that we denote the irreducible parts (i.e. tx^^v^^^) by 
the same letters as in [75j |. but that these quantities here are 
only one half as large as in 75], due to different conventions in 
the definition of torsion. Similarly, the quantities ci,C2,C3 in 
Eq. (1926 are four times as large as in |75H . 



C. Static, spherically and parity symmetric 
vacuum solution 

Hayashi and Shirafuji derived the exact static, spheri- 
cally and parity symmetric i? M „ p(T = vacuum solutions 
for this Lagrangian in [75| . The parallel vector fields take 
the following form in isotropic rectangular coordinates 
(here Latin letters are spatial indices) [751 ]: 



pr 
° = 0, 



toq 
pr 



-p/2 



-l+p/2 



too 
qr 



q/2 



too 
qr 



-1-9/2 



SI , (97) 



where to is a parameter with units of mass and will be 
related to the physical mass of the central gravitating 
body in Section [X] The new parameters p and q are 
functions of a dimensionless parameter e: 

_ k(ci+c 2 ) 



1 + k(c\ + 4c 2 ) 
2 



P 
<1 



1 - 5e 
2 

l-5e 



{{(l~e)(l-4e)]^-2e}, (99) 
{ [ ( l_ e)( l_4 e )]V2 + 2e} . (100) 



Here k = 8ttG. 

The line element in the static, spherically and parity 
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symmetric field takes the exact form [75[ 

V w ) \ or j 

+( i -^)"( 1+ ?) 2+, ^' <101) 

i 



In order to generalize this solution to the axisymmetric 
case, we transform the parallel vector fields into stan- 
dard spherical coordinates and keep terms to first order 
in rrio/r (the subscript "sp" stands for "spherical"): 



3 (sp) k 



I 

k 



(l + ms. 

r 

o 
o 





1 - 


m 


S 1H 






r 


i 




1 - 


mo 


( 1H 






r 







1 - 



sin 8 cos < 
sin 8 sin <; 
cos 8 





COS COS <p 

r 

cos sin <b 



o 

esc sin cp 

r 

esc cos <f) 
r 

o 



(102) 



A particularly interesting solution is that for the pa- 
rameter choice ci = — C2 so that e = and p = q = 2. 
Eq. (jlOip shows that the resultant metric coincides with 
the Schwarzschild metric around an object of mass mo- 
The parameter C3 is irrelevant here because of the static, 
spherically and parity symmetric field. When c\ + C2 is 
small but nonzero, we have e <C 1 and 



P = 2 
q = 2 



e + 0(e 2 ), 
9e + 0(e 2 ) 



(103) 
(104) 



By using equations (f8~4"|) , ([55]) and ([57]) , we find that the 
linearized metric and torsion match our parametrization 
in Section [III Al When e <C 1, the line element is 



ds 2 



dt 2 



1 + 2(1 -2e) 



m 



and the torsion is 

t ™o 

c 9 _ c 4> 
°r9 — D r<t> — 

both to linear order in toq / r. 



"^0 

2r 2 



dr 2 + r 2 dfi 2 , 
(105) 

(106) 
(107) 



D. Solution around Earth 

We now investigate the field generated by a uniformly 
rotating spherical body to first order in e Q . It seems rea- 
sonable to assume that to first order the metric coincides 
with the Kerr-like metric, i.e. 



9t4 = Go(m a a/r) sin 2 



(108) 



around an object of specific angular momentum a in the 
linear regime mo/r <C 1 and a/r <^ 1. Since the Kerr-like 



metric automatically satisfies G({ }) 
vacuum field equation reduces to 



in vacuum, the 



1 



H y.v _ -g^L 2 = . 



(109) 



We now employ our parametrization with "mass" in 
Eq. (|22p replaced by mo, where mo is the parameter in 
accordance with Section IVlII CI In Section[Xl we will ap- 
ply the Kerr solution Q = — 2 after re-scaling mo to corre- 
spond to the physical mass. Imposing the no-curvature 



condition R u 



0, we find that this condition and 



Eq. (|109p are satisfied to lowest order in mo/r and a/r if 



(0) 
(0) 

'A 

(0) 
(0) 

•A 



Go - «o, , 
-2(Q - a ) . 



(o) 



w 



= 2a n . 



a . 



(110) 



Here a superscript (0) indicates the parametrization with 
mo in place of m. ao is an undetermined constant and 
should depend on the Lagrangian parameters c\ , C2 and 
C3. This parameter has no effect on the precession of a 
gyroscope or on any of the other observational constraints 
that we consider, so its value is irrelevant to the present 
paper. 

The parallel vector fields that give the Kerr metric, 
the connection and the torsion (including the spherically 
symmetric part) via equations l[83" l) -(|8i] ) and |86)) - l(87|) 
take the following form to linear order: 

ef = e (sp) k ^ + 



I 
k 



{ r> \ nina sin sin c 

-(£o - a )— — zs - 



(Go - «o) 



mo a sin v cos <; 
^ 



V 



o 



-ao^T\ 




(111 
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IX. A TOY MODEL: LINEAR INTERPOLATION 
IN RIEMANN-CARTAN SPACE BETWEEN GR 
AND HAYASHI-SHIRAFUJI LAGRANGIAN 



We found that the Hayashi-Shirafuji Lagrangian ad- 
mits both the Schwarzschild metric and (at least to linear 
order) the Kerr metric, but in the Weitzenbock spacetime 
where there is no Riemann curvature and all spacetime 
structure is due to torsion. This is therefore an oppo- 
site extreme of GR, which admits these same metrics in 
Riemann spacetime with all curvature and no torsion. 
Both of these solutions can be embedded in Riemann- 
Cartan spacetime, and we will now present a more gen- 
eral two-parameter family of Lagrangians that interpo- 
lates between these two extremes, always allowing the 
Kerr metric and generally explaining the spacetime dis- 
tortion with a combination of curvature and torsion. Af- 
ter the first version of this paper was submitted, Flana- 
gan and Rosenthal showed that the Einstein-Hayashi- 
Shirafuji Lagrangian has serious defects |l26l |. while leav- 
ing open the possibility that there may be other viable 
Lagrangians in the same class (where spinning objects 
generate and feel propagating torsion). This Lagrangian 
should therefore not be viewed as a viable physical model, 
but as a pedagogical toy model admitting both curvature 
and torsion, giving concrete illustrations of the various 
effects and constraints that we discuss. 

This family of theories, which we will term Einstein- 
Hayashi-Shirafuji (EHS) theories, have an action in in 
Riemann- Cart an space of the form 

Ig = J d 4 x^[j-R({}) + cr 2 Cl t x ^tx^ 

+cr 2 c 2 v% + cr 2 c 3 a^a M ] (112) 

where a is a parameter in the range < cr < 1. Here 
the tensors tx^u, and a M are the decomposition (in 
accordance with EqslSBl— f9"0|) of <r~ 1 S u ^\, which is inde- 
pendent of a and depends only on e 1 as per Eq. (|114p . 
The function cr 2 associated with the coefficients ci, C2 
and C3 in Eq. (|112|) may be replaced by any other reg- 
ular function of cr that approaches to zero as a — > 0. 
The metric in the EHS theories is defined in Eq. ([84]) . 
Similar to the Hayashi-Shirafuji theory, the field equa- 
tion for EHS theories is obtained by varying the action 
with respect to the tetrad. The resultant field equation 
is identical to that for the Hayashi-Shirafuji Lagrangian 
(Eq. I93p except for the replacement ci.2,3 — > o 2 c\,2.z- 
Also, the S^ p in Eq. (95|) is replaced by a^S^f. Thus 
the EHS Lagrangian admits the same solution for ej 1 . 
Since the metric is independent of the parameter cr, the 
EHS Lagrangian admits both the spherically symmetric 
metric in Eq. (jlOip and the Kerr- like metric in Eq. (| 1 8 1) , 
at least to the linear order. For the spherically symmet- 
ric metric, the parameter e in Hayashi-Shirafuji theory is 
generalized to a new parameter r in EHS theories, de- 



fined by the replacement c\ t i — > a 2 c\^,'- 

_ Ka 2 (ci+c 2 ) 
7 ' ~ 1 + kct 2 (ci + 4c 2 ) ' 



(113) 



The torsion around Earth is linearly proportional to cr, 
given by the parameter a times the solution in Eq. (|106p 
and pTD]) : 



(114) 



By virtue of Eq. ([7]) (the metric compatibility condition), 
it is straightforward to show that the connection is of the 
form 



T»,. 




ae k p d,e k 



(115) 



EHS theory thus interpolates smoothly between metric 
gravity e.g. GR (cr = 0) and the all-torsion Hayashi- 
Shirafuji theory (cr = 1). If a ^ 1, it is straightforward 
to verify that the curvature calculated by the full connec- 
tion does not vanish. Therefore, the EHS theories live in 
neither Weitzenbock space nor the Riemann space, but 
in the Riemann-Cartan space that admits both torsion 
and curvature. 

It is interesting to note that since the Lagrangian pa- 
rameters c\ and C2 are independent of the torsion pa- 
rameter cr, the effective parameter r is not necessarily 
equal to zero when a = (i.e. , cr 2 ci or cr 2 c2 can be still 
finite). In this case (cr — and yet r / 0), obviously 
this EHS theory is an extension to GR without adding 
torsion. In addition to the extra terms in the Lagrangian 
of Eq. (|112p . the extension is subtle in the symmetry of 
the Lagrangian. In the tetrad formalism of GR, local 
Lorentz transformations are symmetries in the internal 
space of tetrads. Here in this a = 0, r 7^ EHS theory, 
the allowed internal symmetry is global Lorentz transfor- 
mations as in the Weitzenbock spacetime, because tx^, 
Vjj. and a M contain the partial derivatives of tetrads (see 
Ea.mg]). So the cr = and t^O EHS theory is a tetrad 
theory in Riemann spacetime with less gauge freedom. 

Since GR is so far consistent with all known observa- 
tions, it is interesting to explore (as we will below) what 
observational upper limits can be placed on both a and 



X. EXAMPLE: TESTING EINSTEIN 
HAYASHI-SHIRAFUJI THEORIES WITH GPB 
AND OTHER SOLAR SYSTEM EXPERIMENTS 

Above we calculated the observable effects that arbi- 
trary Earth-induced torsion, if present, would have on 
GPB. As a foil against which to test GR, let us now in- 
vestigate the observable effects that would result for the 
explicit Einstein-Hayashi-Shirafuji class of torsion theo- 
ries that we studied in Section IVIII Dl and IIXI 
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Hayashi-Shirafuji with mo 


EHS with mo 


Definitions 


metric 
parameters 


W ' 


-2 


-2 


gtt = -1 - H (0) m /r + O(m /rf 
g rr = l + F^mo/r + O(m /r) 2 




2(1 — 2e) 


2(1 — 2r) 


geodetic 
torsions 


(0) 
*i 


— 1 


— a 


anomalous, S tr * = mo/2r 2 
normal, S r / = S r / = mo/2r 2 


Ao) 


-(l-2e) 


— <t(1 - 2r) 


frame- 
dragging 
torsions 


(0) 

w\ 


So — Qo 


<r(<?0 — «o) 


S r £ = w[ ' (m a/2r ) sin 6 1 
•S'e/ = w< 2^ (moa/2r) sin cos 9 
Sti = «4 0) (m a/2r 2 )sin 2 6> 
St/ = ™i 0) (m a/2r 3 ) sine cos 6» 
= uif moa/2r 4 


(0) 


-2(0o -ao) 


-2a(Q - ao) 


(0) 


ao 


cra 


(0) 


ao 


crao 


(0) 
w 5 


2a 


2crao 



TABLE III: Summary of metric and torsion parameters for General Relativity, Hayashi-Shirafuji gravity and Einstein-Hayashi- 
Shirafuji (EHS) theories. The subscript indicates all parameter values are normalized by an arbitrary constant mo (with the 
units of mass) that is not necessarily the physical mass of the body generating the gravity. The parameter ao in frame-dragging 
torsions is an undetermined constant and should depend on the Hayashi-Shirafuji Lagrangian parameters ci, ci and C3. The 
parameter r, defined in Eq. (|98[) and assumed small, is an indicator of how close the emergent metric is to the Schwarzschild 
metric. The values in the column of Einstein-Hayashi-Shirafuji interpolation are those in the Hayashi-Shirafuji times the 
interpolation parameter a. 





GR 


EHS with autoparallels 


EHS with extremals 


Definitions 


mass 


m 


m = mo 


m = (1 — o)mo 


m — mo 


set by Newtonian limit 


metric 
parameters 


n 


-2 


-2/(1 -a) 


-2 


gtt = -1 - Hm/r + 0{m/r) 2 
g rr = l + Tmjr + 0(m/r) 2 
9t<p = G(ma/r) sin 2 


T 


2 


2(1 - 2t)/(1 - a) 


2(1 - 2r) 


G 


-2 


-2 


-2 


geodetic 
torsions 


ti 





-a /(l -a) 


— a 


anomalous, S tr ' = t\ m/2r 2 
normal, S rg e = S r / = t 2 m/2r 2 


t 2 





- 2t)/(1 - a) 


-cr(l - 2r) 


frame- 
dragging 
torsions 


Wi 





a(G - a) 


a(G - a) 


S r $ = Wi (ma/2r 2 ) sin 2 9 
Se<jy = w 2 (ma/2r) sin 9 cos 9 
S t / =W 3 (ma /2r 2 ) sin 2 9 
S t / = ™4 (ma/2r 3 ) sin#cos6> 
Str* = w 5 ma/2r A 


W2 





-2a{G - a) 


-2a{G - a) 


W3 





aa 


aa 


Wi 





oa 


aa 


W5 





2aa 


2aa 


effective 
torsions 


Hi 





— a 


—a 


fll = (Wl — W2 — W3 + 2W4 + W5)/(—3Q) 
H2 = {Vll —W3+ W5)/ (~G) 







— a 


—a 


bias 


b t 


1 


1 - 4r/3 


1- a -At/2, 


b t = (l+J r + 2t 2 + \v\ti)/S 
b„ = (-a/2)(l + 3/ii-2 M 2) 


K 


1 


(-G/2)(l-a) 


(-0/2)(1-<t) 



TABLE IV: Summary of metric and torsion parameters for Einstein-Hayashi-Shirafuji (EHS) theories of interpolation parameter 
a in autoparallel scheme and in extremal scheme. All parameter values are normalized by the physical mass m of the body 
generating the gravity. The parameter Q and a are related to Qo and ao in Table Hill by Q = (?o/(l — 0) and a = ao/(l — a) in 
autoparallel scheme, G = Go and a = ao in extremal scheme. The value for G is set to —2 by the Kerr metric in linear regime 
m/r <C 1 and a/r <C 1. 





General Relativity 


EHS with autoparallels 


EHS with extremals 


Averaged Geodetic Precession 


(3m/2ro)wo x So 


(1 - 4r/3)(3m/2r )w o x S 


(1 - cr - 4r/3)(3m/2r )w o x S 


Averaged Frame-dragging 


(J/2rg)uJj5 X So 


(-G/2)(l-<j)(I/2rl)u B x So 


(-e/2)(l-ff)(//2rg)aJ B xS 


Second moment (22 


(3/w B /4rg)i x S 


(-3e/w B /8rg)(l - cr)z x S 


(-3e/w B /8rg)(l - cr)i x S - T)(Tmwo(Sj;2 + Sg£)/4r 


Second moment 62 


(3Iue/4t%)x x So 


(-3e/wB/8rg)(l-cr)i x S 


(-3e/w B /8r ( : *)(l - cr)£ x So - ri<rmuo{S%x - Sp)/4r 



TABLE V: Summary of the predicted Fourier moments of the precession rate for General Relativity and the Einstein-Hayashi- 
Shirafuji (EHS) theories in autoparallel scheme and in extremal scheme. 77 = +1 for extremal scheme using S^ v , and —1 for 
extremal scheme using . Other multiple moments vanish. Here m and Iloe are the Earth's mass and rotational angular 
momentum, respectively. 



21 



Effects 


Torsion Biases 


EHS in autoparallel scheme 


EHS in extremal scheme 


PPN biases 


Shapiro time delay 


At/At < - GR 1 = (F-H)/4, 


1 + a-T 


1 - T 


(l + 7)/2 


Deflection of light 


S /s(GR) = (f~H)/4 


1 + cr-T 


1 — T 


(l + 7)/2 


Gravitational redshift 


(Av/v)/{Av/vf GK ) = -H/2 


1 + CT 


1 


1 + a 


Geodetic Precession 


n G /n (GH) = b t 




l-cr- %T 


(1 + 2 7 )/3 


Frame-dragging 


n F /n (GH> = b„ 


1 - a 


1 - cr 


(l + 7 + ai/4)/2 



TABLE VI: Summary of solar system experiments (1): the biases relative to GR predictions for the Einstein-Hayashi-Shirafuji 
(EHS) theories. Both parameters r and a are assumed small. The biases in the PPN formalism are also listed for comparison, 
taken from (3|. 



Effects 


PPN 


EHS in autoparallel scheme 


EHS in extremal scheme 


Remarks 


Shapiro time delay 


7 - 1 = (2.1 ± 2.3) X 10" 5 


cr - r = (1.1 ± 1.2) X lO" 6 


t = (-1.1 ± 1.2) x 10" 5 


Cassini tracking 122] 


Deflection of light 


7 - 1 = (-1.7 ±4.5) x 1CT 4 


a-r= (-0.8 ±2.3) X 10" 4 


T = (0.8 ±2.3) x 10~ 4 


VLBI 123] 


Gravitational redshift 


|et| < 2 x 10" 4 


\a\ < 2 X 10~ 4 


no constraints 


Vessot-Levine rocket |l24j] 


Geodetic Precession 


|7 - 1| < i.i x icr 4 


|r| < 5.7 X 10" 5 


|o- + 4t/3| < 7.6 X 10~ 5 


Gravity Probe B 


Frame-dragging 


7- 1+ \a\\ < 0.024 


\a\ < 0.012 


\a\ < 0.012 


Gravity Probe B 



TABLE VII: Summary of solar system experiments (2): constraints on the PPN and EHS parameters. The constraints on PPN 
parameters are taken from Table 4 and Page 12 of [3|. The full results of Gravity Probe B are yet to be released, so whether 
the frame dragging will agree with the GR prediction is not currently known. The last two rows show the limits that would 
correspond to a GPB result consistent with GR, assuming an angle accuracy of 0.5 milli-arcseconds. 



There are four parameters ci, C2, C3 and cr that define 
an EHS theory via the action in Eq. (|112p . We will test 
EHS theories with GPB and other solar system exper- 
iments. For all these weak field experiments, only two 
EHS parameters — r (defined in Eq. (|113|) ) and cr, both 
assumed small — that are functions of the said four are 
relevant and to be constrained below. 

The predicted EHS metric and torsion parameters, 
studied in Section ITXl are listed in Table Below, 
we will test both the autoparallel and extremal calcula- 
tion schemes. In each scheme, the physical mass m will 
be determined by the Newtonian limit. All metric and 
torsion parameters are converted in accordance with m 
and listed in Table llVl Then the parameter space (r, cr) 
will be constrained by solar system experiments. 



so the physical mass of the central gravitating body is 
m = (1 — cr)mo • (H?) 



Table IIVI lists values of metric and torsion parameters 
in accordance with the physical mass m. Using these pa- 
rameters, the precession rates of gyroscopes in GPB orbit 
can be calculated via equations (|Glj]) . (f6*9"l) . (fT(I)) and ([T5j) . 
The results are listed in Table IVl For GPB, the average 
precession rates are the only experimentally accessible 
observables in practice. GPB will measure the preces- 
sion of gyroscopes with respect to two different axes: the 
orbital angular velocity Coo (geodetic precession) and the 
Earth's rotational angular velocity u>e (frame-dragging). 
As indicated in Table [V] the geodetic precession and 
frame-dragging rates are 



A. Autoparallel scheme 

Hayashi-Shirafuji maximal torsion theory is inconsis- 
tent with the autoparallel scheme, since ti — H/2 = 
(see t\ and Ti in Table ITTT|) . By Eq. ([57| . this means that 
dv/dt = + 0(m/r) 2 . The violation of Newton's law 
rules out the application of the autoparallel scheme to 
the Hayashi-Shirafuji theory. 

However, the Einstein-Hayashi-Shirafuji theories can 
be consistent with this scheme. Using Table Hill the New- 
tonian limit can be written as 



dv 
~dt 



-(1 



m „ 



(116) 



(I - -t)Q 



(GR) 
G 1 



^ = i-|)(i-^r 3 



(118) 
(119) 



where q'q^ and O^?^ are the geodetic precession and 
frame-dragging rate predicted by General Relativity, re- 
spectively. 

The existing solar system experiments, including 
Shapiro time delay, deflection of light, gravitational red- 
shift, advance of Mercury's perihelion, can put con- 
straints on the parameters r and a. The derivation of 
these constraints essentially follow any standard text- 
book of General Relativity [H| except for more general 
allowance of parameter values, so we leave the techni- 
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FIG. 4: Constraints on the EHS parameters (a, r) from so- 
lar system tests in the autoparallel scheme. General Rel- 
ativity corresponds to the black dot (a = t = 0). The 
shaded regions in the parameter space have already been ruled 
out by Mercury's perihelion shift (red/dark grey), the deflec- 
tion of light (orange/grey), Shapiro time delay (yellow/light 
grey) and gravitational redshift (cyan/light grey). If the 
geodetic precession and frame-dragging measured by Grav- 
ity Probe B are consistent with GR to the target accuracy 
of 0.5 milli-arcseconds, this will rule out everything outside 
the hatched region, implying that < a < 8.0 x 10~ 5 and 
-2.3 x 10~ 5 < t < 5.7 x 10~ 5 . Preliminary result of Gravity 
Probe B have only confirmed the geodetic precession to about 
1%, thus bringing no further constraints beyond those from 
gravitational redshift. 



FIG. 5: Constraints on EHS parameters (a, r) from solar sys- 
tem tests in the extremal scheme. General Relativity corre- 
sponds to the black dot (a = r = 0). The shaded regions have 
already been ruled out by Mercury's perihelion shift (red/dark 
grey), the deflection of light (orange/grey) and Shapiro time 
delay (yellow/light grey). If the geodetic precession and 
frame-dragging measured by Gravity Probe B are consistent 
with GR to the target accuracy of 0.5 milli-arcseconds, this 
will rule out everything outside the hatched region, implying 
that < a < f.l x I0~ 4 and -2.3 x 10~ 5 < r < 0.1 x I0" 5 . 
The preliminary result of Gravity Probe B have confirmed the 
geodetic precession only to about 1%, implying that a < 0.01. 



B. Extremal scheme 



cal detail in Appendix [Cl with the results summarized in 
Table ED 

It is customary that biases of GR predictions are ex- 
pressed in terms of PPN parameters on which observa- 
tional constraints can be placed with solar system ex- 
periments. In EHS theories, these biases are expressed 
in terms of the parameters r and a. Thus we can place 
constraints on the EHS parameters r and a by setting up 
the correspondence between PPN and EHS parameters 
via the bias expression. Table PVTl lists the biases in the 
PPN formalism for this purpose, and Table IVlTl lists the 
observational constraints on the EHS parameters r and 
a with the existing solar system tests. 

If GPB would see no evidence of the torsion induced 
precession effects, the (t,ct) parameter space can be fur- 
ther constrained. Together with other solar system ex- 
periments, the observational constraints are listed in Ta- 
ble IVIII and shown in Figure |H 



Einstein-Hayashi-Shirafuji theories predict 7i — —2 re- 
gardless of r and a. By the Newtonian limit, therefore, 
the physical mass of the central gravitating body is just 
the mass parameter mo, i.e. m = tuq. So the parameter 
values do not need rcscaling and are re-listed in Table 
IIVI By these parameters the precession rates can be cal- 
culated and listed in Table |V] As indicated in Table |V] 
the geodetic precession and frame-dragging rates are 

n G = (i-a-^ T )nf B \ (120) 
n F = (i-a)nf R) . (121) 

It is worth noting again that the extremal scheme is not 
a fully consistent framework from the theoretical point 
of view. However, it serves perfectly to show the role of 
EHS theories as the bridge between no-torsion GR and 
Hayashi-Shirafuji maximal torsion theory. Figure [5] illus- 
trates this connectivity in terms of the predictions of GR, 
Hayashi-Shirafuji theory and the intermediate < a < 1 
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FIG. 6: Predictions for the average precession rate by General 
Relativity, Hayashi-Shirafuji (HS) gravity and Einstein-Hayashi- 
Shirafuji theories (for the case of r = and the Kerr solution 
Q = —2) that interpolate between these two extremes, in the ex- 
tremal scheme. Sqp is the geodetic precession rate around the 
orbital angular velocity vector ujo an d #fd is the angular frame- 
dragging rate around Earth's rotation axis lje- The shaded areas 
of about 0.5 milli-arcseconds per year in radius are the approxi- 
mate forecast GPB measurement uncertainties. The two calcula- 
tion schemes using S 1 ^ and 5 M " with extremals for the Hayashi- 
Shirafuji Lagrangian (labeled "HS" in the figure) agree on the 
predicted average rates. The unpublished preliminary results of 
Gravity Probe B have confirmed the geodetic precession to better 
than 1%, so this already rules out the Hayashi-Shirafuji Lagrangian 
and most EHS theories in the extremal scheme in the sense that 
a < 0.01. 



EHS theories, taking r = and Kerr solution Q = —2, 
on the average precession rate (the Sq in Table [Vj . The 
EHS theories are seen to connect the extreme GR and HS 
cases with a straight line. If the data released by GPB 
ends up falling within the shaded area corresponding to 
the GR prediction, the Hayashi-Shirafuji Lagrangian will 
thus have been ruled out with very high significance, and 
the GPB torsion constraints can be quantified as sharp 
upper limits on the cr-parameter. 

More generally, Gravity Probe B will improve the con- 
straints on the (t,ct) parameter space by its precise mea- 
surements of precession rates, in addition to the con- 
straints put by existing solar system experiments. These 
constraints are listed in Table IVlTl and shown in Figure [5] 
As before, the technical details are given in Appendix [C] 



C. Preliminary constraints from GPB's 
unpublished results 

In April 2007, Gravity Probe B team announced that, 
while they continued mining the data for the ultimately 
optimal accuracy, the geodetic precession was found to 
agree with GR at the 1% level. The frame-dragging yet 
awaits to be confirmed. Albeit preliminary, these unpub- 
lished results, together with solar system tests, already 



place the first constraint on some torsion parameters to 
the 1% level. More quantitatively, \t 2 + < 0.01 

in the model-independent framework, while W\ + w% — 
w 3 — 2w4 + w 5 is not constrained. In the context of 
EHS theories, the constraint is scheme dependent. In the 
autoparallel scheme, GPB's preliminary results place no 
better constraints than those from gravitational redshift 
(~ 10~ 4 ). In the extremal scheme, however, the prelim- 
inary results give the constraint a < 0.01. The bottom 
line is that GPB has constrained torsion parameters to 
the 1% level now and will probably reach the 10 -4 level 
in the future. 



XI. CONCLUSIONS AND OUTLOOK 

The PPN formalism has demonstrated that a great way 
to test GR is to embed it in a broader parametrized class 
of theories, and to constrain the corresponding parame- 
ters observationally. In this spirit, we have explored ob- 
servational constraints on generalizations of GR includ- 
ing torsion. 

Using symmetry arguments, we showed that to lowest 
order, the torsion field around a uniformly rotating spher- 
ical mass such as Earth is determined by merely seven di- 
mensionless parameters. We worked out the predictions 
for these seven torsion parameters for a two-parameter 
Einstein-Hayashi-Shirafuji generalization of GR which 
includes as special cases both standard no-torsion GR 
(a = 0) and the no-curvature, all torsion (er = 1) 
Weitzenbock spacetime. We showed that classical so- 
lar system tests rule out a large class of these models, 
and that Gravity Probe B (GPB) can further improve 
the constraints. GPB is useful here because this class 
of theories suggested that, depending on the Lagrangian, 
rotating objects can generate torsion observable with gy- 
roscopes. In other words, despite some claims in the 
literature to the contrary, the question of whether there 
is observable torsion in the solar system is one which ul- 
timately can and should be tested experimentally. 

Our results motivate further theoretical and experi- 
mental work. On the theoretical side, it would be inter- 
esting to address in more detail the question of which 
Lagrangians make torsion couple to rotating objects. A 
well-defined path forward would be to g eneralize the 
matched asymptotic expansion method of jllll Ill2| to 
match two generalized EHS Kcrr-likc Solutions in the 
weak-field limit to obtain the laws of motion for two well- 
separated rotating objects, and determine which of the 
three non-equivalent prescriptions above, if any, is cor- 
rect. It would also be interesting to look for generaliza- 
tions of the EHS Lagrangian that populate a large frac- 
tion of the seven torsion degrees of freedom that symme- 
try allows. Finally, additional observational constraints 
can be investigated involving, e.g., binary pulsars, grav- 
itational waves and cosmology. 

On the experimental side, Gravity Probe B has now 
successfully completed its data taking phase. We have 
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shown that the GPB data constitute a potential gold 
mine of information about torsion, but that its utility for 
constraining torsion theories will depend crucially on how 
the data are analyzed and released. At a minimum, the 
average geodetic and frame dragging precessions can be 
compared with the predictions shown in Figure [6j How- 
ever, if it is technically feasible for the GPB team to 
extract and publish also different linear combinations of 
the instantaneous precessions corresponding to the sec- 
ond moments of these precessions, this would enable look- 
ing for further novel effects that GR predicts should be 
absent. In summary, although the nominal goal of GPB 
is to look for an effect that virtually everybody expects 
will be present (frame dragging), it also has the potential 
to either discover torsion or to build further confidence 
in GR by placing stringent limits on torsion theories. 
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APPENDIX A: PARAMETRIZATION OF 
TORSION IN THE STATIC, SPHERICALLY AND 
PARITY SYMMETRIC CASE 



In this appendix, we derive a parametrization of the 
most general static, spherically and parity symmetric tor- 
sion in isotropic rectangular and spherical coordinates. 
The symmetry conditions arc described in Section [III A II 
with the quantity O now being the torsion tensor S„ v p . 
Note that torsion (the antisymmetric part of the connec- 
tion) is a tensor under general coordinate transformations 
even though the full connection is not. 

First note that time translation invariance is equivalent 
to the independence of torsion on time. Then consider 
time reversal, under which a component of torsion flips 
its sign once for every temporal index. Invariance un- 
der time reversal therefore requires that non-zero torsion 
components have either zero or two temporal indices. To- 
gether with the fact that torsion is antisymmetric in its 
first two indices, this restricts the non-zero components 
of torsion to be S 0i ° and (i = 1, 2, 3). 

Now consider the symmetry under (proper or im- 
proper) rotation (see Eq. (TT3|) ). The orthogonality of 
the matrix R enables one to write 



dx' 1 
~dxJ 



R ij , 



dx 1 



at' 



at 



Thus formal functional invariance means that 

& Oi ( X ') = R %J &0j ( X ) = ( X ')i mm 

S^ix') = W m R kn R a S mn l {x) = S jh ?(x<). 



Eq. (|A2|) requires that the torsion should be built up 
of x l and quantities invariant under 0(3), such as 
scalar functions of radius and Kronecker ^-functions, 
since 8[, r = /»"''/>'•''•' = R 1 ' 1 ^' 1 = R i ' i {Rr y )^' = 
6i>ji . Note that we are interested in the parity sym- 
metric case, whereas the Levi-Civita symbol e^jf. is a 
three-dimensional pseudo-tensor under orthogonal trans- 
formations, where "pseudo" means that is a ten- 
sor under SO(3) but not under 0(3), since (! itl y k t = 

R l l Ri j R k k eijk = deti? x evjiy ■ Therefore, eyfe is pro- 
hibited from entering into the construction of the torsion 
tensor by Eq. (|A"2|) . 

Thus using arbitrary combinations of scalar functions 
of radius, x l and Kronecker ^-functions, the most general 
torsion tensor that can be constructed takes the form 



c o 

°0i 



, 771 i 

2r 6 



J 3i) i 



(A3) 
(A4) 



where the combinations t\m and are arbitrary func- 
tions of radius. Note that in Eq. (|A4|) . terms proportional 



to x l x J x or x l 5jk are forbidden by the antisymmetry of 
the torsion. We will simply treat the functions ti(r) and 
t?{r) as constants, since GPB orbits at a fixed radius. 

Transforming this result to spherical coordinates, we 
obtain 



c t 



<? - 



t dx 1 m 

j dxi dx k d9 m 
° k dr 86 dx^^ h 2^ 



; dx? dx k 



00 m 
—— = fa — 

dx 1 2r 2 



at 



, 1. (Al) 



All other components not related by the antisymmetry 
vanish. In the above equations, the second equalities 
follow from the chain rule and the facts that dx 1 /dr = 
x % = e£, dx 1 /d9 — reg, and dx 1 /d<j> = rsin^e^, where e*, 
e l s and are the zth-components of the unit vectors in 
spherical coordinates. To first order in the mass m of the 
central object, we need not distinguish between isotropic 
and standard spherical coordinates. 



APPENDIX B: PARAMETRIZATION IN 
STATIONARY AND SPHERICALLY 
AXISYMMETRIC CASE 

Above we considered the Oth order contribution to the 
metric and torsion corresponding to the static, spheri- 
cally and parity symmetric case of a non-rotating spheri- 
cal source. In this appendix, we derive a parametrization 
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of the most general 1st order correction (denoted by a 
superscript (1)) to this metric and torsion that could be 
caused by rotation of the source, i.e. corresponding to the 
stationary and spherically axisymmetric case. The sym- 
metry conditions are described in Section IIII B 11 with 
the quantity O replaced by the metric gl^J for Appendix 
IB II and by the torsion S$ p for Appendix IB 21 



In spherical polar coordinates where the z-axis is par- 
allel to J, this first-order correction to the metric takes 
the form 

g2?=0^*?8, (B5) 

where ma — J is the magnitude of J. All other compo- 
nents vanish. 



The Metric 



The Torsion 



The invariance under time translation makes the met- 
ric time independent. Under time reversal J — * —J, and 
a component of the metric flips its sign once for every 
temporal index. Thus, the formal functional invariance 
equation for time reversal reads 



±g$(x\J)=g$(x\-J). 



(Bl) 



The plus sign in Eq. (|B1[) is for components with 
even numbers of temporal indices, and minus sign for 
those with odd numbers. Since only terms linear in 
J/r 2 = e m e a are concerned, the minus sign in the argu- 
ment —J can be taken out as an overall factor, implying 
that the non-vanishing components of metric can have 
only one temporal index. Thus the only nonzero first- 
order correction to in rectangular coordinates is 
(i=l,2,3). 

Now consider the transformation property under 
(proper or improper) rotation. By the orthogonality of 
the matrix R, the vector x transforms as x — > x' = Rx 
(Eq. (jAip V Since J is invariant under parity, formally 
the transformation of J writes as 



J -> J' = (<MR) x RJ . 



(B2) 



The formal functional invariance for rotation reads 



&g§\x\J) 



9«' & 



V) 



(B3) 



That J is a pseudo-vector under improper rotation re- 
quires that the Levi-Civita symbol e^-fc, also a pseudo- 
tensor, appear once and only once (because J appears 
only once) in the metric so as to compensate the detH 
factor incurred by transformation of J. Other possible el- 
ements for construction of the metric include scalar func- 
tions of radius, x l , J\ &y. Having known the elements, 
the only possible construction is therefore 



fit 



(i) 



(B4) 



where x l = x l jr is the unit vector of position vector and 
Q is dimensionless. Assuming that there is no new scale 
other than the angular momentum J built into the 1st 
order of torsion theory, i.e. no new dimensional param- 
eter with units of length, Q(r) must be a constant by 
dimensional analysis, since the factor J 1 has explicitly 
appeared. 



We follow the same methodology as for our 
parametrization of the metric above. Given the time- 
independence, the property that J reverses under time- 
reversal requires that the non-vanishing components 
of torsion have only one temporal index, so they are 
S'^j*, S^lj (i,j=l,2,3) in rectangular coordinates. (The 
antisymmetry of torsion over its first two indices excludes 
the possibility of three temporal indices.) Under (proper 
or improper) rotation, the formal functional invariance 
equation reads 



s ( %VlJ) 



R^R^S {1 k Y(x\3) = S ( gV|J'), 



Sg?V|J) = R*B? l S$>(x\J) = S£]W). 

Again, in building the torsion, one should use the Levi- 
Civita symbol e^fc once and only once to cancel the detR, 
factor from the transformation of J. The most general 
construction using scalar function of radius, x" 1 , 5ij, J 1 
(also appearing once and only once) and is 

? (i)t _ J]_ 
' ij ~ 2r 3 ' 



= —< llk J k + ^J k x l (e,u.r 



r-3 



£jklX l ) , 



J:: "' 1 A J k x l e lkl x 3 + I^J k x l e lk ix l 



2r 3-v- 1 2r 3 2r 3 

By the same dimensional argument as in Appendix (|B 1[) , 
fx, . . . , /s must be dimensionless constants. 

Transforming the above equations to spherical coordi- 
nates where the z-axis is parallel to J, we obtain to first 
order 

, dx l dx 3 ma . , „ 



^ °« dr d(j) 

c (l) t Q t 0x1 dxj 



13 39 



ma . 
W2 sin cos ( 



2r 



o(l) r 
° t<fr 

c-(l) 6 



dx % dx 3 
d(j) dr 



gg dx 1 dx 3 ma 



2r 2 



S 



(i) <t> 



9 Stij-z T ~5 7T = -rsin0COS0, 
)k oO 2r 6 

dx 1 dx 3 



dr dcf> 



y tij ~d6 d(j> 

All other components vanish. The constants are related 
by Wi = fx - f 2 , w 2 = ft, w 3 = U ~ w 4 = -/ 3 , 
w 5 = h + h- 
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Planet 



Earth 




Sun ^0 

FIG. 7: Geometry of the Shapiro time delay measurement. 



APPENDIX C: CONSTRAINING TORSION 
WITH SOLAR SYSTEM EXPERIMENTS 



1. Shapiro time delay 

For the electromagnetic field, if torsion is coupled to 
the vector potential by the "natural" extension, i.e. 
, d^Ay — > W fj,A v using the full connection, the Maxwell 
Lagrangian — jF^F'- 11 ' will contain a quadratic term in 
A^ that makes the photon massive and breaks gauge 
invariance in the conventional form. Since the photon 
mass has been experimentally constrained to be < 10~ 17 
eV, we assume that A^ does not couple to torsion. In- 
stead, we assume that the Maxwell field Lagrangian in 
the curved spacetime with torsion follows the extension 
d^Ay — > vj? A v using the Levi-Civita connection. Since 
the Levi-Civita connection depends on the metric and its 
derivatives only, light rays follow extremal curves (metric 
geodesies). 

In general, assume the line element in the field around 
a (physical) mass to is 



tional field of the Sun is 





[■r=D 


/•v—De 


T = 2 


\ dt + 


/ dt 




Jr=D p 


Jr=D 



2[Jm - d 2 



D% - D 2 ] + H)i 



x In 



-J-m 



UD* -D 2 + D p )(y/D%-D* + D E ) 



D 2 




(C3) 



If D < D E and D < D p , the third term in Eq. (|C3|) is 
negligible compared to the second one. The excess travel 
time At of a round-trip light ray is 



At = T-2[JD$-D 2 



where At^ GR ^ is the excess time predicted by GR 



At^^ = 4m In 



(De + x E ■ n)(D p — Xp-n) 



(C4) 



(C5) 



Here x E (x p ) is the vector from the Sun to the Earth 
(the planet), and n is the unit vector from the planet to 
Earth (see Figure [7j). 

For EHS theories in the autoparallel scheme, [T — 
Ti.)/4 = (1 - e)/(l - a) w 1 + a - e, if a < 1. For 
EHS theories in the extremal scheme, (T — 7i)/4 = 1 — e. 



ds 2 



1 



H— 

r 



dt 2 



1 



T— 

r . 



dr 2 +r 2 dn 2 . (CI) 



The effect of the rotation of the mass can be ignored 
when the rotation is slow. 

Light deflection angle is tiny for the solar system tests 
we consider, so a ray can be well approximated by a 
straight line. Let us use coordinates where the Sun (of 
mass to), the Earth and a planet reflecting the light ray 
are all in the x-y plane (8 = n/2) and the cc-axis points 
along the ray from the planet to Earth (see Figure [7]). 
Let D be the minimal distance of the ray from the Sun. 
Then r sin0 = D, or rd(f> = — tan^rfr. Since ds 2 = for 
a light ray, 



dt 2 



dt 



(1 



w !^)-i(i + ^!!! + tan 2 ) dr 2 

r r 
lr 2 r , to mD 2 , 

— [i + (^-w)7-^-pH 



r\dr\ , to rnD 2 . 



(C2) 



The round-trip travel time for an electromagnetic signal 
bouncing between Earth and the Planet in the gravita- 



2. Deflection of light 



As discussed in Appendix IC 1[ we assume that a light 
ray follows an extremal curve (metric geodesic), taking 
the form 



Dt 




dx v dx p 
dr dr 



= 0. 



(C6) 



Here /Dt denotes the covariant differentiation using 
the Levi-Civita connection. 

The fx = t component of the metric geodesic is 

d 2 t dt dr ^ 

dr 2 r 2 dr dr 

or, to order 0(m/r), where to is the mass of the Sun 
deflecting the light, 



d_ 



771, dt 

(l + H-) — 

r dr 



0. 



Integrating this gives a conserved quantity, 

const . 



TYl , dt 

k=(l + H—)— 
r dr 



(C7) 
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The \i = 9 component of the metric geodesic admits the 
planar solution 9 = tt/2. The fi = (f> component of the 
metric geodesic, when 9 = 7r/2, is 

d 2 4> 2 dr dej) 



dr 2 



r dr dr ' 



whose first integral gives another conserved quantity, 



— „2 ( 



dr 



const . 



(C8) 



For light rays in the equatorial plane 9 = 7r/2, 



ds^ 
dr 2 



1 + H- 



dr 
(C9) 

Note that the fi = r component of the metric geodesic is 
not independent of Eq. (|C9[) . Rewriting dt/dr and dtfi/dr 
in terms of k and h via Eq. (|C20|) and Eq. (|C2ip , respec- 
tively, and using the fact that dr/dr — (dr / d(f) (d<j) / dr) , 
one finds 



4fi 



u = —Trau 
2 



k 2 T + TL 



-m , 



h 2 2 

where u = 1/r. The solution to order 0(m) is 



(CIO) 



sin< 



Tm 



(1 + Ccos0 + cos 2 



k 2 T + n 



D 2D 2 T YJ h 2 2 

(Cll) 

where D is the minimal distance of the ray to the Sun. 
The a;-axis is set up to be along the incoming direction 
of the ray. C is an arbitrary constant that can be deter- 
mined at 4> = 7r (incoming infinity) . As long as deflection 
angle (5 <g; 1, 



2Tm k 2 _ „„„ 
— -J^ m (F + H)D. 



Using 



k 



2#/, 



w-) 

r 



(C12) 



(C13) 



is the angular momentum of the light ray relative to the 
Sun, we finally obtain 



S (GR) 



(C14) 



where S^ GR ^ = Am/D is the deflection angle predicted by 
GR to lowest order. 



3. Gravitational Redshift 

As discussed above, we assume that the orbits of light 
rays are metric geodesies even when there is non-zero tor- 
sion. Non-relativistically, the metric geodesic equation 
for a test particle is 



dv 
~dt 



(-H)m 



(C15) 



Effectively this introduces the gravitational potential U, 
defined by dv/dt = F= -VZ7, 



U 



(-H) m 



2 r 

Thus the gravitational redshift of photons is 

(GR) 



(C16) 



(C17) 



where (Ai/ / 'u)^ GR ^ is the redshift predicted by GR 



v 



(GR) 



m , 1 1 



(C18) 



For EHS theories in the autoparallel scheme, —Ti/2 = 
1/(1 — cr) m 1 + crforcr^; 1. For EHS theories in extremal 
scheme, — 7Y/2 = 1 exactly. 



4. Advance of Mercury's Perihelion in autoparallel 
scheme 



In the autoparallel scheme, a massive test particle (e.g. 
a planet in the field of the Sun) follows an autoparallel 
curve (i.e. an affine geodesic). We now derive the advance 
of the perihelion when torsion is present. The autopar- 
allel equation reads 



Du» d 2 x t " 



Dt 



dr 2 



dx^dx^ = 
"P dr dr 



(C19) 



where D/Dt is the covariant differentiation by the full 
connection. 

The /i = t component of Eq. ()C19|) reads 



dr 2 



(ti-ny 



n dt dr 
■ 2 dr dr 



= 0, 



or, to order 0(m/r), where m is the mass of the central 
gravitating body (e.g. the Sun), 



d_ 



I* ,„ , .m\ dt 
(l + (H-h)j) 



dr 



0. 



The integral gives a conserved quantity k, 

m\ dt 
r / dr 



k= M + (ft -th- 



eorist . 



(C20) 



The \i = 9 component of Eq. (|C19[) admits the planar 
solution 9 = tt/2. The fi = (f> component of Eq. (|C19I) . 
when 9 = tt/2, is 



d 2 (j> 2 ^ m dr d(f> 
dr 2 r r 2 dr dr 



0, 



whose first integral gives another conserved quantity h, 

2 d( t ) l 



V--T- l+t 2 - 

dr r 



const . 



(C21) 
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The path parameter r can be chosen so that 

dx>* dx v 



ds /dr = g^ v - 



dr dr 



-1 . 



(C22) 



Eq. (1C22|) is consistent with the autoparallcl scheme since 
V p g M „ = and Du^/Dt = 0. Note that the fi — r com- 
ponent of Eq. (|C19[) is not independent of Eq. (IC22|1 . For 
a test particle in the equatorial plane 9 = tt/2, Eq. (|C22|) 
reads 



„ .m 
l + H— 

r 



l + T— 

r 



.dr. 2 7,d<j).n 



1 . 



(C23) 

Reusing the trick employed in Appendix IC 2| we find 



(C24) 



d U 3^. o HI r . 9 , ^ 

— 2+ u= -Tmu^-^-U-T- 
+2h + 2t 2 ) + F - 2t 2 ] , 



to order 0(mu), where u = 1/r. Note that to lowest 
order k « 1 + C(to, (velocity) 2 ) , so the second term on 
the right hand side of Eq. (|C24j) becomes (t 1 -H/2)m/h 2 . 
Since m is the physical mass of the central gravitating 
body, the autoparallel scheme requires t\ — H/2 = 1. 
Now Eq. (|C24|) becomes 



d 2 u 



m 3 , 
- + -Tmu . 



(C25) 



Solve the equation perturbatively in the order of e 



(m/h) , i.e. use the ansatz u 



u 

Ul 



h 2 

2h 2 



(1 + ecos< 



wo + eu\. One finds 

(C26) 
(C27) 



e z 1 
1 + ecj) sin <f) + — (1 — - cos 2</>) 



Eq. (|C26|) gives the classical elliptical orbit with eccen- 
tricity e and the semi-latus rectum p = a(l — e 2 ) = h 2 /m. 
The sin </> term in Eq. (|C27[) contributes to the advance 
of the perihelion, while the constant and cos 2<f> terms do 
not. Therefore 



m 
h 2 



ecos 



•iTra 2 
~2h 2 ~ 



(C28) 



In Eq. (|C28[) , we used the fact that the second term inside 
the cosine is <C 1. The advance of the perihelion is now 
given by 



AO = 



2tt 



1 - 
T 



2h? 



2tt 



(C29) 



where A8^ GR ^ = 6irm 2 /h 2 — (titm/p is the perihelion 
advance predicted by GR. 



5. Advance of Mercury's Perihelion in extremal 
scheme 

The extremal scheme assumes that a test particle (e.g., 
a planet) follows the metric geodesic even though the 
torsion is present. Following the same algebra as in Ap- 
pendix IC 41 and noting that Ti = — 2 for the extremal 
scheme, we finds that the advance of the perihelion in 
the extremal scheme has the same bias factor J 7 / 2, i.e. , 
Eq. (IC291) holds. 



APPENDIX D: CONSTRAINING TORSION 
PARAMETERS WITH THE UPPER BOUNDS 
ON THE PHOTON MASS 

In this Appendix, we derive the contraints on torsion 
parameters that result from assuming that the "natural" 
extension — > V M (using the full connection) in the elec- 
tromagnetic Lagrangian. This breaks gauge invariance, 
and the photon generically gains a mass via an additional 
term of the form —^m^g^A^A^ in the Lagrangian as we 
will now show. The assumption gives 

/„„ - 2S^ X A X , (Dl) 

d v A^. The Maxwell Lagrangian 



where = d^A. 
therefore becomes 



C 



EM 



!J y r p.v r afi , 



g^g^Uvf*? - K^ApA,, + S^ x A x f^ , 

(D2) 



where K^ v = S a p S a > . The Euler-Lagrange equation 
for the action S — J d^x^Z—gC^M yields the following 
equation of motion for A^: 



V}J^ = 2S^f Xu + 2K Xu A x + 2Vj}(S^ x A x ) . (D3) 

Here Vj^ and vj^ are the covariant derivative w.r.t. the 
full connection and the Levi-Civita connection, respec- 
tively. Both the 2nd and 3rd terms on the right hand 
side of Eq. (|D3[) contain the coupling to A^. To clarify 
this, Eq. (|D3j) can be rewritten non-covariantly as 

v£/"" = 2S fi /f Xl/ + 2A\ [K Xu + d^S^ 




2S^ x d„A> 



(D4) 



in which the 2nd term on the right hand side is the direct 
coupling of A^ . 

The matrix is symmetric. If it is also positive 
definitive up to the metric signature ( — h H — h), the first 
term in the square bracket may be identified as the pho- 
ton mass term. In the field of a non-rotating mass, using 
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the parametrization (Eqs. [T4l and [T5|) . it can be shown 
that 



2r 4 ' 



(D5) 
(D6) 
(D7) 



The matrix K has the eigenvalues — 2pri ^ (with eigen- 

vector f) and | 4 (with 2 degenerate eigenvectors). 
Since the metric signature is ( — h++), all photon masses 
are positive or zero, The nonzero ones are of order 



m 

m 7 ~ t— 



or (with units reinserted) 



2 fiG m 

m^c ~ t 7t . 

c r z 



(D8) 



(D9) 



Here t = max(\tx\, fa\) and r is the distance of the ex- 
periment location to the center of the mass m that gen- 
erates the torsion. For a ground-based experiment here 
on Earth, this gives 



t ~ 4.64 x 10 22 m 7 c7(l eV) . 



(D10) 



The upper bound on the photon mass from ground-based 



experiments is m 7 c < 10 eV [1251 ] . so the constraint 



that this bound places on the dimensionless torsion pa- 
rameters is quite weak. 

Experimentalists can also search for an anomalous elec- 
tromagnetic force and translate the null results into pho- 
ton mass bounds. To leading order, the anomalous force 
is 2d^S^ uX A\, since the K-term is proportional to S 2 , 
while the 2nd term in the square bracket of Eq. (|D4|) is 
proportional to S. In a field of a non-rotating mass m, 



uvX\00 



(d^s^y* = t 



(c^S*^) * = {d„S^ x Y° = 0,(D11) 



2r 3 



x l x 3 



(D12) 



which has eigenvalues x (0,-1,-1,2). This cannot 
be identified as a mass term since there must be a nega- 
tive "mass squared" regardless of the sign of t%. However, 
the anomalous electromagnetic force expressed as a pho- 
ton mass can be estimated as 



fa\h 2 G- 



or 



\t 2 \ a 1.23 x 10 18 m 7 c 2 /eV. 



(D13) 



(D14) 



This implies that current ground-based experimental up- 
per bounds on the photon mass are too weak (giving 
merely \t\ < 10 2 , as compared to \t\ = 1 from Hayashi- 
Shirafuji gravity) to place constraints on torsion param- 
eters that are competitive with those from GPB. 
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